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AN  F-LAYER  MODEL  OF  THE  IONOSPHERE  FOR  NORTH  AMERICA 


I 


by 

C.J.  Gladwin 


ABSTRACT 


<r 


Data  collected  from  the  Alouette  1  and  2  satellites 
(ALOSYN  and  ALOUETTE  N(h>)  have  been  organized  in 
order  to  make  predictions  of  ionospheric  character¬ 
istics  such  as  foF2,  slab  thickness,  and  spread  F 
as  functions  of  longitude,  latitude,  local  time, 
month,  magnetic  activity,  and  sunspot  number. 

Spread  F  an<l_f^F2  are  determined  statistically 
from  about&lO^^louette-l  ionograms  ajjd  slab 
thickness  is  estimated  using  about  ^O^/electron 
density  profiles.  The  data  and  the  necessary 
computer  programs  for  analysing  these  data  by 
the  user  are  on  two  labelled  tapes  and  are 
available  from  the  Communications  Research  Centre. 
Contour  maps  of  foF2,  spread  F  and  slab  thickness 
as  a  function  of  geographic  latitude  and  local 
time  can  be  generated.  A  special  technique  is 
outlined  in  order  that  the  user  may  modify  his 
results  to  include  the  main  trough.  The  major 
computer  programs  used  are  available  on 
request . 


o 


A 


1.  INTRODUCTION 

The  Alouette  1  Topside  Sounder  satellite  was  launched  on  September  29, 
1962.  From  that  date  until  end  of  1968,  a  continuing  effort  of  the  Defence 
Research  Board  Telecommunications  Establishment  (DRTE,  the  forerunner  of  the 
Communications  Research  Centre,  CRC)  was  to  routinely  scale  the  approximately 
1,000,000  ionograms  recorded,  for  several  ionospheric  characteristics  such  as 
the  F2  layer  penetration  frequency,  foF2,  as  well  as  a  subjective  spread 
index  in  the  neighbourhood  of  that  point.  The  data  have  been  stored  on 
magnetic  tapes  and  have  come  to  be  known  as  ALOUETTE  1  IONOSPHERIC  DATA, 
ALOSYN  [1], 
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Concurrently,  60,000  of  these  lonograms  as  well  as  another  20,000 
Alouette  2  lonograms  were  scaled  for  "almost-complete"  topside  electron 
density  profiles  (that  is  within  about  3%  of  the  foF2).  These  data  likewise 
have  been  stored  on  magnetic  tapes  and  become  known  as  ALOUETTE  1  (or  2) 
IONOSPHERIC  DATA,  N(h)  [2], 

The  intent  of  this  report  is  to  describe  a  procedure  for  organizing  and 
reducing  these  vast  quantities  of  data  in  order  that  a  user  can  obtain  predic¬ 
tions  of  various  ionospheric  characteristics  (foF2,  spread,  slab  thickness, 
and  possibly  MUF(3000)F2)  as  functions  of  longitude,  latitude,  local  time, 
sunspot  number,  month,  and  magnetic  activity  (as  measured  by  the  3-hour ly  Kp 
index).  It  should  be  noted  that  while  MUF(3000)F2  was  not  a  scaled  parameter, 
it  can  be  computed  from  an  estimated  height  of  the  maximum  electron  density 
by  means  of  an  empirical  formula  given  by  Shimazaki  [27]  or  Schmid  et  al  [26], 
In  addition  to  the  prediction  of  various  characteristics,  the  model  has  the 
capabilities  of  computing  percentiles  for  the  characteristics.  A  technique 
is  also  presented  to  enable  the  user  to  estimate  regions  of  severe  density 
gradients  such  as  the  main  trough  [17],  It  is  intended  that  the  model  also 
be  used  to  supplement  the  work  of  Muldrew  and  James  [19]  in  which  the 
ionospheric  effect  on  the  Doppler  frequency  of  a  Search  and  Rescue  satellite 
was  estimated. 

Global  ionospheric  modelling  using  vast  quantities  of  data  came  of  age 
in  the  early  1960's  with  the  advent  of  easy  accessibility  to  large  computers. 
The  first  major  one  was  the  C.C.I.R.  model  which  is  accessible  either  through 
an  atlas  of  ionospheric  characteristics  [5]  or  a  companion  punched  card 
version.  On  a  smaller  scale  (geographic  latitude  >  35°N) ,  there  exists  the 
Canadian  model  [23]  and  that  of  Petrie  [20],  Essentially  these  models  are 
based  on  monthly  medians  of  scaled  characteristics  for  each  hour  of  local 
time;  the  latter  being  observed  at  about  100  ground  stations  around  the 
globe.  The  time  frame  was  the  five  years  in  the  mid  1950's  corresponding  to 
the  rising  part  of  the  sunspot  cycle.  Using  the  numerical  mapping  procedures 
developed  by  Jones  and  Gallet  [13],  these  data  were  reduced  to  sets  of  monthly 
coefficients  for  each  characteristic,  and  then  these  sets  were  normalized  to 
two  sunspot  numbers  (Rz**0  and  Rz*»100) .  The  predictions  were  assumed  to  vary 
linearly  with  sunspot  number  for  Rz  in  [0,  150]  and  be  independent  of  it 
otherwise. 

The  next  extensive  modelling  exercise  was  culminated  by  Llewellyn  and 
Bent  [15].  Again  they  predicted  foF2  from  a  set  of  bottomside  stations 
(400,000  lonograms),  but  in  addition  they  used  the  60,000  Alouette  1 
profiles  and  6,000  electron  density  probe  results  (1967-68)  to  construct  a 
"closed  form"  density  profile.  The  inputs  for  their  model  were  longitude, 
latitude,  time,  season,  solar  flux,  sunspot  number,  and  foF2.  Their  closed 
form  profile  essentially  ignored  underlying  layers  and  ledges  (FI,  E,  and  D) 
because  they  were  primarily  interested  in  obtaining  accurate  values  for  total 
electron  content,  TEC,  only.  Thus  their  bottomside  TEC  could  be  used  in 
conjunction  with  the  present  model,  in  place  of  the  approximation  to  be 
outlined  later. 

For  discussion  and  listings  of  some  other  modelling  experiences  one 
should  see  Kohnlein  [14],  Another  model,  which  will  be  applied  here,  is  one 
of  the  trough  created  by  Halcrow  and  Nisbet  [10], 
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The  model  to  be  outlined  In  this  report  will  now  be  considered.  While 
there  are  1,000,000  topside  ionograms  the  vast  majority  of  these  is  limited 
to  the  continents  of  North  and  South  America;  see  Figure  1.  However,  unlike 
the  data  base  for  the  bottomslde  models,  there  exists  virtually  continuous 
latitudinal  coverage.  Due  to  the  precession  of  the  orbital  plane  of  Alouette 
1  and  the  revolution  of  the  earth  about  the  sun,  observations  at  a  given 
latitude  are  obtained  for  two  4  hour  local  time  intervals  per  month.  The 
local  time  history  of  Alouette  1,  coincidentally,  repeats  itself  almost 
exactly  every  six  months.  Hence,  the  local  time  is  essentially  the  same  at  a 
given  geographic  latitude  for  a  given  day  of  any  year.  Thus,  in  making 
predictions  for  a  given  month,  at  least  3  months  of  data  will  be  required. 

On  the  other  hand  instead  of  using  monthly  medians  of  characteristics,  semi¬ 
monthly  medians  will  be  used  because  of  the  nature  of  the  local  time  cells. 

A  few  other  minor  nuisances,  arising  essentially  from  attempting  to  design 
an  empirical  model  after  the  data  have  been  collected,  are  also  noted: 

(1)  Because  the  satellite's  solar  cells  degraded,  the  number  of 

observations  for  the  last  3  years  of  the  recording  process  was 
virtually  halved  (see  Table  2).  Thus,  there  will  be  a  "bias"  in 
the  model  towards  the  lower  sunspot  years  1963-66. 


Figure  1.  Location  and  Coverage  of  Alouette  1  Telemetry  Stations.  After  [<?] 
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(2)  Because  of  the  routine  scaling  procedures,  and  possible  change  in 
habits  of  the  scalers  over  the  6  year  period,  the  number  of 
observations  (of  foF2  or  fxf2)  which  are  "unreliable'1  (i.e.  scaled 
to  an  uncertainty  of  more  than  0.2  MHz)  can  range  up  to  60%  in  any 
given  month  (see  Table  2) .  This  problem  of  routine  scaling  is  not 
as  serious  for  the  N(h)  data,  however,  as  these  ionograms  were 
scaled  on  user  demand. 

(3)  No  data  are  available  over  the  geographic  poles  (due  to  the  80° 
inclination  of  the  orbital  plane  of  Alouette  1). 

(4)  The  spread  index  may  not  be  as  consistent  as  that  obtained  in 
smaller  scale  spread  experiments  [4], 

(5)  The  maximum  sounding  frequency  of  Alouette  1  is  11.5  MHz.  Hence, 
this  model  cannot  give  accurate  foF2  values  if  the  expected  value 
exceeds  about  10  MHz.  The  model  will  often  give  foF2  values 
greater  than  11.5  MHz  due  to  extrapolation  to  high  sunspot  numbers. 

Only  a  model  for  North  America,  where  most  of  the  data  were  measured, 
is  considered  (see  Figure  1).  The  actual  domain  for  the  model  will  be: 

(1)  geographic  longitude  in  [-180°,  -45®],  that  is  [180°W,  45°W], 

(2)  geographic  latitude  in  (-10®,  80°],  that  is  [10°S,  80°N],  or 
geomagnetic  latitude  approximately  in  [0®,  90®], 

(3)  all  sunspot  numbers,  using  CCIR  assumptions  for  Rz  >  150, 

(4)  all  local  times,  months  and  Kp. 

Then  for  a  given  month  and  year,  there  will  be  3  equidistant  longitude 
zones  and  18  equidistant  geomagnetic  latitude  zones,  as  well  as  4  zones  of 
local  time  and  Kp,  each.  Initially  the  year  is  kept  separate  because  of  the 
varying  sunspot  number. 

An  interesting  comparison  of  results  is  provided  by  Petrie  and 
Lockwood  [21]  where  they  used  geographic  longitude  zones  of  30  degrees, 
geographic  latitude  zones  of  2  degrees,  and  local  time  zones  of  1  hour  for 
ALOSYN  data  from  1963. 

The  computer  used  in  this  study  is  a  Honeywell  Sigma  9  operating  under 
Control  Program  -  Five  (CP-V) .  The  data  analysis  programs  and  user  programs, 
available  on  request,  are  written  in  Fortran  IV  and  hence  should  be  adaptable 
to  most  computers. 

In  chapter  2,  a  description  of  the  ALOSYN  data  base  is  given  and  the 
implementation  of  the  prediction  model,  for  foF2  and  spread,  is  considered. 

A  similar  discussion  is  taken  up  in  chapter  3  for  the  ALOUETTE  N(h)  data 
base  in  order  to  compute  quantities  such  as  TEC  and  slab  thickness.  Chapter 
4  is  concerned  with  Integrating  the  trough  model  [10]  with  the  model  described 
here.  The  notations  used  are  summarized  in  Appendix  B.  A  glossary  of  files 
and  tapes  is  given  in  Appendix  A  relating  the  various  programs  to  their 
respective  sections  in  the  report.  A  brief  description  of  the  files  referred 
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to  in  the  text  is  also  given  in  Appendix  A.  A  source  listing  of  the  files  is 
also  available  on  request. 


2.  ALOSYN  MODEL 


2.1  DESCRIPTION  OF  THE  DATA  BASE 


The  entire  collection  of  1,000,000  ionograms  is  stored  on  8  9T  (9  track) 
tapes  in  year-universal  time  order.  They  have  been  blocked  into  3200  byte 
records.  Each  record  contains  at  most  40  ionograms  with  80  bytes  of  informa¬ 
tion.  A  new  record  is  started  with  each  day  change.  There  is  an  80  byte 
header,  or  title,  to  be  skipped  at  the  beginning  of  each  tape.  A  header  has 
also  been  placed  at  the  beginning  of  some  months.  A  complete  description  of 
all  the  data  in  a  record  may  be  found  in  Nelms  et  al.  [3]  or  at  the  beginning 
of  any  of  the  monthly  or  semi-monthly  ALOSYN  booklets  [1]  that  have  been 
published.  Here  only  a  list  of  the  data  actually  used  in  the  model  is 
presented: 


Year 

Month 

Day 

Universal  time  (hour,  minute,  second) 
Geographic  local  time  (hour,  minute) 
Geographic  longitude  (degrees) 

Geographic  latitude  (degrees) 

Geomagnetic  latitude  (degrees) 

Height  of  the  satellite  (kilometres) 

Solar  angle  at  the  satellite  (degrees) 
Gyrofrequency  at  the  satellite  (MHz) 
Penetration  frequency  of  the  0  wave  (MHz) 
Accuracy  of  scaling  the  0  wave 
Spread  of  the  0  wave 

Penetration  frequency  of  the  X  wave  (MHz) 
Accuracy  of  scaling  the  X  wave 
Spread  of  the  X  wave 
Planetary  magnetic  activity  index 


Y 

M 

D 

Ut 

Lt 

Long 

Lattg 

Lattm 

hs 

*s 

w 

f  oF2 
a(o) 
s(o) 
fxF2 
a(x) 
s(x) 

Kp 


In  addition,  one  should  consult  Table  3  which  displays  the  smoothed 
monthly  Zurich  sunspot  numbers  which  have  been  used.  Most  of  the  above  items 
are  self-explanatory;  also  see  the  accuracy  and  spread  definitions  which  may 
be  found  in  Tables  1(a)  and  1(b)  respectively. 
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TABLE  1(a) 

Accuracy  Convention  of  ALOSYN  foF2 


a  (o) 


1 

2 

3 

4 

5 

-0 


foF2 

foF2 

foF2 

foF2 

foF2 

foF2 


Definition 


within  .05  MHz  of  tabulated  value 
within  .1  MHz  of  tabulated  value 
within  .2  MHz  of  tabulated  value 
less  than  tabulated  value 
greater  than  tabulated  value 
not  observed 


TABLE  1(b) 

Spread  Convention  of  ALOSYN 


S 


A, 

D, 

0, 

J. 


(o) 


Definition 


Index 


B,  C 
E,  F 
H,  I 
K,  L 


No  spread 
Slightly  spread 
Moderately  spread 
Extremely  spread 


1 

2 

3 

4 


In  Table  1(a)  there  is  some  doubt  about  the  interpretation  of  a«4  or  5, 
due  to  changes  in  scaling  habits  over  the  years.  In  the  early  years  these 
appear  to  have  meant  an  uncertainty  greater  than  0.2  MHz  and  a  -0  was  used 
for  "unobservable"  critical  frequencies.  However,  the  correspondence  given 
in  the  table,  a  in  [1,5],  is  adopted  for  all  years.  Using  this  table,  foF2 
is  defined  for  each  ionogram,  as  follows  (see  Appendix  B) . 


foF2 


(2.1) 


[fxF2(fxF2-f  (hmax))  J*5  ,  for  a(x)  In  [1,3] 
n 

foF2,  for  (a(o),a(x))  in  [1,3]  x  [4,5] 

where  f„(hmax)  is  found  by  an  inverse  cube  extrapolation  from  hg  to  hmax  = 

300  km. 

In  Table  1(b)  the  three  letters  associated  with  each  group  were  used  to 
record  other  information  (e.g,  oblique  traces,  cusps,  etc.),  but  the  index  at 
the  right  will  be  sufficient  here.  In  a  similar  vein,  Kp  values  defined  as 
n-,  no,  or  n+  are  taken  to  be  n  -  .3,  n,  or  n  +  .3  where  n  is  an  integer  in 
[0,9]. 

TABLE  2 

Blocked  ALOSYN  Records  and  Percentage  of  Unreliable  lonograms  by  Month  and  Year 


"\VYEAR 

month\. 

1962 

1963 

1964 

1965 

1966 

1967 

1968 

JANUARY 

507 

485 

458 

2  28 

235 

136 

10 

28 

45 

42 

29 

44 

FEBRUARY 

__ 

mmm 

439 

412 

301 

230 

166 

mm ■ 

26 

52 

56 

26 

42 

MARCH 

_ 

503 

428 

463 

198 

268 

171 

28 

25 

48 

30 

26 

41 

APRIL 

562 

554 

571 

MM 

285 

213 

19 

18 

38 

■El 

20 

26 

MAY 

621 

mm 

MM 

334 

269 

214 

13 

■El 

in 

18 

20 

34 

JUNE 

493 

418 

502 

328 

285 

215 

6 

21 

21 

18 

21 

34 

_ 

592 

582 

■Mgm 

261 

271 

228 

10 

13 

mmm 

12 

_ _ _ 

18 

22 

AUGUST 

_ 

615 

495 

590 

■ 

266 

203 

11 

26 

23 

■E9 

19 

32 

SEPTEMBER 

20 

392 

433 

230 

256 

172 

16 

WEM 

31 

24 

22 

24 

31 

OCTOBER 

613 

602 

316 

561 

258 

274 

164 

34 

28 

32 

28 

20 

25 

38 

NOVEMBER 

637 

599 

406 

190 

199 

■■ 

29 

26 

■El 

38 

32 

37 

m 

DECEMBER 

506 

479 

419 

271 

173 

143 

28 

26 

62 

46 

27 

50 

wEmm 

TABLE  3 

(Smoothed)  Sunspot  Numbers,  Rz 


w 


jiiU.iW&i 
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In  Table  2,  the  following  statistics  are  displayed: 

(a)  the  actual  number  of  (blocked)  records  with  data,  by  month  and 
year, 

(b)  the  percentage  of  "unreliable"  ionograms,  by  month  and  year, 

i.e.  those  for  which  an  foF2  is  not  definable,  using  equation  (2.1). 
For  example,  for  May  1965,  there  are  600  blocked  records  and  18%  of 
the  total  number  of  ionograms  are  unreliable. 

In  view  of  the  introductory  remarks,  the  actual  number  of  ionograms  in 
a  given  month  and  year  is  approximately 

I  =  40  R  -  20  Dm  (2.2) 

where  R  =  number  of  blocked  records 

Dm  =  number  of  days  in  month  M 

I  =  actual  number  of  ionograms. 


2.2  DEVELOPMENT  OF  THE  MODEL 

2.2.1  Organization  of  the  Data  Base 

A  description  of  how  the  data  were  processed  is  given  here  and  in 
Figure  2.  A  description  of  how  the  user  uses  these  processed  data  is  given 
in  section  2.3. 


Figure  2.  Schematic  Overview  of  Model  Development 


Instead  of  processing  the  8  tapes  sequentially  (i.e.  year  by  year),  they 
were  processed  month  by  month  for  all  years.  While  this  increases  the 
requisite  number  of  tape  usages  by  (at  most)  a  factor  of  76/8  =9.5,  it 
permits  a  modular  processing  procedure  for  the  vast  quantities  of  data  (see 
below).  Also,  this  procedure  is  more  natural  in  the  sense  that  convenient 
checks  arise  on  the  consistency  of  the  data.  Thus,  the  first  step  is  to  set 
up  Nyear  files  for  a  given  month,  ALOMONXX,  where  Nyear  is  the  number  of  years 
of  data  available  for  that  month,  and  XX  are  the  last  two  digits  of  the 
respective  year  represented  by  ALOMONXX.  Each  of  these  files  contains  R  +  1 
records  (see  Table  2),  i.e.  each  file  contains  either  the  header  for  the 
month  or  the  last  record  of  the  previous  month,  followed  by  the  R  records  in 
the  month.  During  this  initial  processing  stage,  in  addition  to  completing 
Table  2,  counts  of  the  number  of  complete  or  partial  records  omitted  due  to 
"tape  problems"  are  also  kept.  This  is  necessary  in  order  to  position  the 
tape  correctly,  for  starting  the  next  month.  Additional  checks  are  provided 
by  actually  printing  out  the  first  4  records  as  well  as  the  last  3  of  each  of 
the  files  ALOMONXX.  The  control  files  AXX  and  BXX  are  then  executed  to 
continue  the  next  stage  of  processing  (see  below). 

2.2.2  Geographic  Smoothing  and  Interpolation 

In  the  first  stage  of  this  process,  executing  AXX,  the  ionograms  are 
sorted  into  a  number  of  longitudinal  and  latitudinal  zones  of  essentially 
constant  magnetic  activity  and  local  time.  The  longitude  range  is  divided 
into  3  equal  zones  and  the  geomagnetic  latitude  range  is  divided  into  18 
equal  zones  (see  the  model  domain  given  in  chapter  1) .  Each  of  these  54 
geographic  areas  (see  Figure  3)  is  then  subdivided  into  4  magnetic  activity 
zones  defined  by  the  ranges  of  Kp: 

[0,1.3]  ;  [1.7, 3. 3]  ;  [3. 7, 5.3]  ;  [5. 7, 9.0]; 

and  4  local  time  zones  (see  below  for  a  more  complete  description) .  Thus 
there  are  864  cells,  although  many  of  these  will  be  empty  in  the  southwest 
corner  (over  oceans,  see  Figure  3,  or  Table  4).  Furthermore  the  fourth  Kp 
zone  (corresponding  to  severe  storms)  is  fairly  scant  in  terms  of  data.  How¬ 
ever,  the  fine  division  is  used  to  prevent  features  such  as  the  main  trough 
from  becoming  too  "smoothed  out";  the  latter  is  strongly  dependent  on  local 
time  and  Kp  [17]. 

These  zones  define  864  cells.  Associated  with  each  of  the  864  cells, 
a  record  is  kept  consisting  of  the  following  data  array: 

(n^,  Long^,  Lattm^,  x^,  Kp^»  Lt^,  foF2max^,  foF2min^)  k  =  1(1)864, 

where  n.  is  the  number  of  ionograms  in  the  cell  and  the  other  numbers  are 
arithmetical  averages,  minima  or  maxima  of  the  respective  observations  in 
the  cell.  These  data  are  placed  on  a  scratch  file  DUMXX  which  is  updated  on 
the  next  pass  through  ALOMONXX  by  BXX.  This  update  consists  of  appending  a 
frequency  distribution  consisting  of  10  equidistant  intervals  of  [foF2minjt, 
foF2maxjt]  for  foF2,  and  also  a  distribution  for  the  spread  (with  4  intervals). 
The  choice  of  10  intervals  for  foF2  was  partially  motivated  by  the  number  of 
ionograms  n^  being  tabulated  and  also  by  the  fact  that  the  range  (*  foFlmax^  - 
foF2min  )  is  roughly  2  MHz,  so  that  the  grouping  process  has  an  uncertainty 


of  about  the  same  magnitude  as  the  maximum  scaling  uncertainty  (corresponding 
to  accuracy  a  =  3) . 

This  updated  collection  of  864  records  is  put  on  a  file  called  MONXX 
together  with  the  Long,  Lattm,  Lt,  Kp,  M,  and  Y  subscripts  (see  below).  This 
completes  the  first  reduction  which  has  seen  an  original  ALOSYN  month  reduced 
by  a  factor  of  about  20  in  computer  storage.  In  the  next  section  a  further 
reduction  by  a  factor  of  about  10  (by  eliminating  empty  cells  and  standardizing 
to  one  sunspot  number)  will  be  carried  out. 

A  brief  description  of  the  construction  of  the  4  local  time  zones  and 
their  respective  indices  is  now  given. 


Figure  3.  Geographic  Zone  Alignment  of  the  ALOSYN  Mode I 
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In  Figure  4,  three  orbits  of  the  satellite  have  been  plotted  (geographic 
latitude  v.  geographic  local  time)  for  the  1st,  16th  and  31st  of  October.  For 
other  months,  see  the  caption.  The  4  zones  are  described  by  a  pair  of  sub¬ 
scripts  (LT1,LT2)  defined  on  the  product  set  [1,2]  x  [1,2],  as  shown  in  Figure 
4.  The  first  subscript  defines  the  "half-month"  that  the  particular  ionogram 
is  in  and  is  determined  by  comparing  the  actual  hour  in  the  month  (Greenwich) 
of  the  ionogram  to  the  total  number  of  hours  in  the  month.  The  second  sub¬ 
script  determines  whether  the  ionogram  is  from  a  northbound  (daytime,  in  this 
case)  pass  or  not.  This  is  determined  by  whether  or  not  the  local  time  of  the 
ionogram  is  less  than  the  theoretical  local  time  at  the  maximum  geographic 
latitude  of  the  satellite  (approximately  80°N),  Ltmax.  The  latter  has  been 
computed  assuming  a  linear  variation  with  year  fraction,  Yf  (where  Yf,  for 
example,  for  January  31  1963,  is  about  63.085). 

Ltmax  =  16.48  -  47.93(Yf  -  62.74)  ±  .05  hours,  modulo  24.  (2.3) 

This  empirical  formula  was  determined  by  passing  a  straight  line  through 
the  actual  (Ltmax,  Yf)  pairs  from  1962  day  272  at  2130  hours  and  1966  day  274 
at  0120  hours.  At  these  times  Alouette  1  was  actually  near  80°N  and  the 
actual  values  of  Ltmax  were  computed  by  fitting  a  parabola  to  5  (Lattg,  Lt) 
pairs  in  the  neighbourhood  of  80°N. 


Figure  4.  Definition  of  the  Local  Time  Zones  for  October;  April  is  the  same.  For  November  and  May, 
December  and  June,  January  and  July,  February  and  August  and  March  and  September,  Move  the 
Curves  4,8, 12, 16 JO  hours,  respectively,  to  the  Left. 
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Actually,  the  local  time  zones  are  not  so  symmetrically  placed  about 
Ltmax  since  geomagnetic  latitude  is  used  in  defining  the  Lattm  zones.  Thus 
the  4  local  time  zones  deviate  from  the  expected  2  hour  zone  size  (especially 
around  the  poles),  but  this  was  considered  a  reasonable  "trade-off"  in  view 
of  the  well  known  alignment  of  foF2  contours  with  magnetic  latitude  (e.g.  see 
[17]). 


Finally,  for  calculating  the  average  local  time  of  cell  k  (to  ensure 
continuity  around  midnight),  the  following  recursive  formula  is  used: 

A±  =  [ (i-1)  (A±_x  +  24r)  +  (Lt1  +  24s)]  /  i  (2.4a) 

where  i  =  1(1)0^,  r  and  s  are  integers  in  [0,1]  such  that  the  quantity 
| (A^  +  24r)  -  (Lt^  +  24s) |  is  a  minimum.  The  average  local  time  for  cell  k 

is 


Lt,  -  A  (2.4b) 

k 

A  cell  (or  record)  of  any  MONXX  shall  be  referenced  according  to  seven 
subscripts  LO,  LA,  LT1,  LT2,  LK,  LM,  LY)  defined  on  the  (extended)  product 
set 

[1,3]  x  [1,18]  x  [1,2]  x  [1,2]  x  [1,4]  x  [1,12]  x  [2,8] 

corresponding  to  the  Long,  Lattm,  Lt  (two  values),  Kp,  M,  and  Y  zones, 
respectively.  The  Long  zones  are  numbered  from  west  to  east,  the  Lattm  zones 
from  north  to  south,  the  Lt  zones  as  in  Figure  4,  the  Kp  zones  with  increas¬ 
ing  magnetic  activity,  and  the  M  zones  from  January  to  December.  The  Y  zone 
is  simply  the  second  digit  of  "XX"  in  MONXX.  After  the  standardization  of 
the  MONXX  (see  below),  the  year  subscript,  Y,  will  be  dropped. 

A  single  example  of  a  record  of  MONXX  is  now  given.  A  cell  was  chosen 
where  the  range  of  the  foF2  distribution  was  considerably  greater  than  2  MHz, 
in  order  to  illustrate  the  geographic  smoothing  involved  in  the  model. 

Example  l:  Cell  1,5, 2, 1,3, 2, 7  ;  February  1967. 

The  contents  of  the  record  are  shown  in  Figure  5. 

See  also  Example  2  to  follow. 

On  actually  examining  the  7  ionograms  which  comprised  the  cell  it  was 
noticed  that  the  4  ionograms  at  the  lower  end  of  the  foF2  range  came  from  a 
locale  (165°W,  68°N)  at  1730  hours,  while  the  upper  3  came  a  day  later  and 
an  hour  earlier  at  the  locale  (136°W,  66°N).  Thus  the  two  extreme  points  have 
been  replaced  by  a  hypothetical  intermediate  point,  yielding  a  median  foF2  of 
4.5  MHz. 

An  extract  from  MONXX  starting  at  cell  (3, 7, 2, 1,1, 9, 2)  is  shown  in 
COMPUTER  OUTPUT  1. 
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Figure  5.  Distribution  of  foF2.  for  a  Given  Set  of  Conditions. 


2.2.3  Sunspot  Standardization 

Here  a  procedure  for  reducing  the  Nyear  files  MONXX  to  a  single  file 
MON,  of  essentially  the  same  form,  but  standardized  to  some  central  sunspot 
number  depending  on  the  respective  sunspot  numbers  of  the  individual  MONXX, 
will  be  described.  The  reason  for  selecting  sunspot  number  (rather  than  any 
of  the  other  5  independent  variables)  is  that  ionospheric  characteristics, 
e.g.  foF2,  have  been  observed  to  possess  strong  linear  relationships  with 
sunspot  number,  at  least  up  to  Rz  =  150.  See  Davies  [6],  for  example. 

There  are  a  number  of  ways  in  which  this  standardization  procedure  may 
be  achieved;  here  the  approach  adopted  in  the  C.C.I.R.  model  [5]  is  briefly 
described  to  note  some  of  the  difficulties  in  applying  it  to  satellite 
observations.  This  provides  some  motivation  for  the  procedure  which  was 
actually  adopted. 

For  each  of  the  60  months  (1954-58),  the  data  for  each  hour  of  local 
time  at  approximately  100  ground  stations  have  been  used  to  determine  a 
numerical  map  of  the  form: 


6 

foF2  -  a  +  E  a„.  , 
0  J-l  2j'1 


cos  jT  +  a2j 


sin  jT 


(2.5) 


where  T  =  15  Lt  -  180  is  the  local  time  relative  to  noon,  in  degrees. 
Actually,  each  of  the  a.  are  dependent  on  76  coefficients  to  describe  the 
longitudinal  and  latitudinal  variations  so  that  for  each  month  (of  60)  there 
are  988  (  ■  76  x  13  )  coefficients  in  all.  Then  for  a  given  month,  (all 
years)  and  for  each  of  the  988  coefficients,  a  least  squares  line  is  fitted 
to  the  coefficient  versus  sunspot  number.  From  these  relations,  two  sets  of 


COMPUTER  OUTPUT  1.  Extract  from  a  File  MONXX.  XX  =  62(1)68. 
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988  coefficients  are  tabulated,  for  each  month,  corresponding  to  the  sunspot 
numbers  Rz“0  and  Rz=100. 

For  a  topside  model,  72  months  of  data  (1963-68)  would  have  to  be  fitted 
to  a  map  of  the  above  form  in  3  month  stages  at  least  (to  ensure  complete 
local  time  coverage).  This  would  mean  that  the  aj  would  also  have  to  depend 
on  an  extra  variable,  month,  and  perhaps  even  a  Kp  variable.  This  would  give 
24  (seasonal)  sets  of  approximately  988  coefficients.  Then  for  a  given 
season  and  coefficient,  a  least  squares  line  would  be  fitted  to  the  coeffici¬ 
ent  versus  "seasonal"  sunspot  number  relationship.  There  can  be  quite  a  bit 
of  variation  of  sunspot  number  with  season  —  especially  during  the  years 
1966-68  where  the  data  are  scantiest  (see  Table  2). 

The  procedure  adopted  here  will  be  an  "inverse"  of  the  above  -  i.e. 
predictions  are  made  for  foF2  for  a  given  set  of  conditions  (Long,  Lattm,  Lt, 
Kp,  and  M  from  a  given  record  of  MONXX)  and  then  these  (at  most)  Nyear  numbers 
(one  from  each  year)  are  standardized  to  an  intermediate  sunspot  number. 

The  Nyear  files  MONXX  are  first  merged  into  one,  MONB96,  which  has  a 
blocked  record  format.  That  is,  each  record  consists  of  the  same  cell  for 
all  years  (864  blocked  records).  Unfortunately  the  number  of  ionograms  varies 
considerably  from  year  to  year  for  a  given  cell.  Probably  an  even  more 
serious  difficulty  is  that  in  trying  to  keep  only  "quality"  cells  (with  Nyear 
or  Nyear-1  years  of  data)  the  geographic  coverage  is  drastically  reduced 
overall.  This  is  illustrated  in  Table  5  for  2  particular  months. 

For  brevity,  any  such  cell  will  be  denoted  as  type  (M,N).  The  M  years 
are  not  necessarily  the  first  M  years  (starting  from  1962  or  1963),  but  in 
most  cases  this  will  probably  be  true.  Thus,  for  example,  if  only  cells  with 
at  least  2  years  containing  at  least  2  ionograms  are  desired,  the  "practical 
minimum",  only  447  cells  out  of  595  (75%)  are  kept  for  February. 

The  procedure  actually  adopted  is  to  use  only  the  (2,2)  cells  of 
MONB96.  Thus  for  each  of  these  cells  a  least  squares  fit  is  performed  to  the 
Nyear  pairs  of  median  foF2  and  actual  sunspot  number  to  obtain: 

f oF2  =  a  Rz  +  b  (2.6) 

Then,  in  later  interpolation  and  prediction  (see  section  2.3),  only  those 
cells  which  yield  positive  values  of  a  and  b  are  used.  It  is  noted  here  that 
of  the  10,368  (  ■  12  x  864)  cells  there  are  5,618  of  type  (2,2).  Of  these, 

586  gave  non-positive  values  for  a  or  b.  Thus,  for  use  in  interpolation  and 
prediction,  (see  section  2.3)  there  will  be  5,032  (or  approximately  half  of 
the  total  possible).  As  a  matter  of  interest,  it  is  noted  that  there  are 
1,005  cells  (all  months)  which  only  contain  1  year  of  data.  No  attempt  has 
been  made  to  Incorporate  these  cells  into  the  model  since  coefficients  a  and 
b  in  (2.6)  cannot  be  determined,  independently,  for  each  cell.  Furthermore 
there  are  126  cells  (all  months)  which  yielded  values  of  a  >  0.1.  The 
majority  of  these  are  of  "dubious"  quality  arising  from  calculating  a  and  b 
with  2  or  3  years  at  essentially  the  same  sunspot  number. 

While  some  of  these  coefficients  (a  and  b)  may  be  "rough",  leading  to 
"rough"  values  of  foF2,  these  values  will  be  smoothed  out  later  (see  section 
2.3).  This  "roughness"  and  the  effect  of  smoothing  are  illustrated  in 


A 


TABLE  5 

,N).  The  Number  of  Blocked  Records  (Cells),  for  February  and  October,  which  have  at  Least  M  Years  with  at  Least  N  lonograms  each  Year. 
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Figures  6(a)  and  6(b).  The  local  time  of  these  plots  is  about  1730  hours 
at  4S°N  geomagnetic  latitude.  The  actual  regression  model  for  the  smoothed 
gradients,  a  (MHz/unit  sunspot  number),  is: 

a  =  0.053  -  0.011  sin(irLt/12)  -  0.015  cos(TrLt/12) 

+  0.016  cos2 (TrLt/12)  -  0.037  (Lattm/90)  (2.7) 

+  0.015  (Lattm/90)4  -  0.046  (Kp/9.3)3, 


with  a  root  mean  square  error  of  +  0.019  in  a. 

The  equation  (2.7)  was  determined  by  applying  a  stepwise  regression 
procedure  [24]  to  all  February  cells  which  yielded  values  of  a  and  b.  Only 
statistically  significant  (non-zero)  coefficients  (up  to  order  4)  are 
displayed  in  (2.7). 

The  actual  details  of  the  standardization  of  a  given  cell  of  MONB96, 
assuming  there  are  a  "suitable"  number  of  years  containing  a  "suitable" 
number  of  ionograms  are  now  given.  One  should  also  note  Figure  7  where  the 
mathematical  assumptions  underlying  the  standardization  procedure  are 
illustrated,  namely: 

(1)  the  gradient  of  foF2  with  sunspot  number,  a,  is  positive, 

(2)  the  intercept  at  Rz  =  0,  b,  is  positive, 

(3)  the  conditional  probability  distribution  of  foF2  for  a  given  Rz 
is  independent  of  Rz. 


The  following  notations  are  used  in  the  standardization  procedure.  Let 
(n,  Longk,  Lattmk,  Ltk,  Kpk,  Rzk,  Fining,  Fmaxk,  (A..  :  i  =  1(1)10,  Bki  :  i  = 
1(1)4))  k  =  1(1)  Nyear  2  denote  a  particular  blocked  record  (cell)  of 
M0NB96.  Here  Ak^  and  Bk^  ar^  the  histogram  occurrence  frequencies  on  [Fining, 
Fmaxk]  and  [0.5,  4.5]  for  foF2  and  spread,  respectively.  Note  that 
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First  the  medians  M^(F)  and  Mk(S)  are  computed  from  the  foF2  and  spread 
histograms  using  inverse  linear  interpolation  between  the  right-most  abscissae 
and  the  cumulative  histogram  probabilities.  Using  the  pairs  (Rzk,  Mk(F))  k  ■ 
l(l)Nyear,  regression  coefficients  a  and  b  such  that  foF2  =  aRz  +  b  are 
computed.  Then  using  Rz  =  Ny®ar  Rzk/Nyear  the  standardized  limits  of  the  foF2 

histograms  are  computed: 

SFmin  =  min 
k 


(Fmin^  +  a(Rz  -  Rz^)) 


and  SFmax  =  max 

k 

The  Aki  observations  Fki  =  Rnink  +  (i  -  0.5)hk  where  hk  =  (Ftaaxk  -  Fmink)/10 
are  standardized  to  SFi  =  Ffcl  +  a(Rz  -  Rzk)  and  inserted  into  a  histogram  with 


(Fmax^  +  a(Rz  -  Rz^)) 


Figure  6(b).  Behaviour  of  foF2-sunspot  Number  Gradient,  1.5  <Kp<  3.5. 
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Figure  7.  Assumptions  Concerning  the  foF2  Distribution  (for  a  Given  Cell,  all  Years)  Used  in  the  Sunspot 

Standardization  Procedure. 


10  equidistant  intervals  on  [SFmin,  SFmax]  with  occurrence  frequencies  SA^. 


Note  that 


10  Nyear 

l  SA  =  I  n.  •  N. 
i-1  1  k-1  k 


In  a  similar  fashion,  one  could  standardize  the  spread  distributions. 
However,  due  to  the  subjectiveness  of  that  index  they  have  been  simply 
accumulated  for  all  years.  Essentially  then,  it  is  assumed  that  spread  is 
independent  of  sunspot  number.  As  a  matter  of  interest.  Figure  8(a)  shows 
median  spread  plotted  against  sunspot  number  for  3  different  cells.  The 
local  time  of  these  cells  is  approximately  0900  hours  in  January  and  July  and 
1300  hours  in  December. 

The  spread  occurrence  frequencies  are  then: 

Nyear 

SB,  =  Z  B  ,  i  =  1(1)4 
1  k=l 

Finally  arithmetical  means  of  Long,  Lattm,  Lt,  and  Kp,  for  all  years 
are  computed  for  each  cell. 


A  "keyed”  record  is  written  on  the  file  MONSTD  to  represent  the  stand¬ 
ardized  cell.  The  contents  of  this  record  include:  (Key,  N,  Nyear,  Long, 
Lattm,  Kp,  Lt,  Rz,  SFmin,  SFmax,  SA^,  SBj ,  a,  b),  i  =  1(1)10,  j  =  1(1)4.  (See 
also  COMPUTER  OUTPUT  4,  section  2.3). 

The  keys  associated  with  each  record  are  of  a  type  associated  with 
direct  access  READ/WRITE  statements,  which  enable  the  user  to  directly  access 
standardized  cells  upon  specification  of  his  input  conditions.  The  keys  are 
integers  defined  as  a  linear  combination  of  the  cell  subscripts,  with  the 
weights  depending  on  the  particular  cell  alignment  used.  In  this  case: 

4  LK 
4 

or  more  compactly. 

Key  =  864  LM  +  288  LO  +  16  LA  +  8  LT1  +  4LT2  +  LK  -  1180  (2.8) 

Thus  the  keys  are  in  the  range  [1,  10368]  with  empty  cells  (no  data)  corres¬ 
ponding  to  non-existent  keys.  The  last  step  in  preparing  the  reduced  data 
file  for  the  model  is  simply  to  concatenate  all  the  months  together  using 
PCL  (file  management  system  for  the  Honeywell  computer)  on  a  file  ALLMON. 

A  single  example  is  now  presented. 

Example  2:  February:  1,5, 2, 1,3, 2;  this  is  a  (3,5)  cell  (see  Table  5). 

Here  there  are  3  quality  years*.  1964,  1966  and  1967. 

The  blocked  record  of  M0NB96  appears  as: 
k  nk  Longk  Lattn^  X.tfc  KPr  Fmi,^  FmaXj^  Rzk  A^  BR1 

1  18  -157.0  67.3  16.94  4.0  2.07  5.05  16.2  1480112001  06120 

2  7  -153.0  67.6  17.13  4.6  4.53  5.36  32.2  3200100001  01  60 

3  7  -153.4  67.7  17.08  4.2  2.48  8.36  75.7  2101000111  04  30 

In  Figure  8(b),  the  median  foF2  v.  Rz  relationship  for  these  three 
years  has  been  plotted  for  this  cell.  The  standardized  histogram  of  foF2 
lies  in  [1.71,  7.59]  with  occurrence  frequencies  consisting  of  the  array 

2,  2,  12,  2,  3,  6,  2,  1,  1,  1 

from  minimum  foF2  to  maximum  foF2. 

As  a  consistency  check  it  is  noted  that  for  the  immediately  preceding 
cell  (lower  Kp),  the  foF2  v.  Rz  relationship  is  even  "flatter": 

foF2  -  0.019  Rz  +  3.64. 

This  cell  is  of  much  better  quality,  (5,15),  in  terms  of  Table  5. 


Key  = 


10368 
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MEDIAN  SPREAD  INDEX  IN  [l,4] 
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SUNSPOT  NUMBER ,  Rz 

Figure  8(a).  Spread  Versus  Sunspot  Number  for  a  Given  Longitude,  Latitude,  Local  Time,  Magnetic 

Activity  and  Month. 

GEOGRAPHIC  LONGITUDE  154.  5°W 

GEOMAGNETIC  LATITUDE  67  5°N 

LOCAL  TIME  1703  hrs 

MAGNETIC  ACTIVITY,  Kp  4.3 

FEBRUARY 

NUMBER  OF  I0N0GRAMS  <3yrs.)  32 


0  50  100 

SUNSPOT  NUMBER,  Rz 


Figure  8(b).  foF2  Versus  Sunspot  Number  for  a  Given  Longitude  Latitude,  Local  Time,  Magnetic 

Activity  and  Month. 


f 


24 


2.3  IMPLEMENTATION  OF  THE  MODEL 

2.3.1  Real  Time  Accession  of  the  Data 

The  reduced  file,  ALLMON,  now  exists  which  consists  of  at  most  10,368 
keyed  records,  each  record  representing  a  standardized  cell.  Actually  about 
one  third  of  these  records  will  be  "non-existent"  due  to  empty  cells  over 
oceans  areas. 

Upon  specification  by  the  user  of  a  particular  point  (Long,  Lattg,  Rz, 
M,  Kp),  the  subscripts  of  that  cell  can  be  computed,  the  key  computed  by 
equation  (2.8),  and  the  cell  directly  accessed.  Actually  only  the  shortened 
equation 


Key  =  864  LM  +  288  LO  +  16  LA  +  8(2)  +  4(2)  +  4  -  1180  -  16 

=  864  LM  +  288  LO  +  16  LA  -  1168  (2.9) 

need  be  used  and  (at  most)  the  next  16  cells  for  each  of  the  months  M-l,  M, 
and  M+l.  As  mentioned  above,  3  months  are  required  to  obtain  sufficient 
local  time  coverage.  See  Table  6  for  a  display  of  the  cell  (key)  structure 
of  ALLMON.  With  little  additional  effort,  one  can  access  5  months  (i.e.  also 
M-2  and  M+2);  however,  unless  stated  explicitly  otherwise,  it  will  be  assumed 
throughout  the  remainder  of  the  chapter  that  3  months  have  been  accessed. 

A  few  comments  concerning  the  use  of  Table  6  are  appropriate: 

(1)  This  represents  the  "format"  for  the  month  of  January  (other 
months  are  simply  translations  of  the  key  by  864) . 

(2)  The  3  longitude  zones  are  represented  by  the  3  major  columns,  and 
the  18  latitude  zones  by  the  18  major  rows. 

(3)  Within  one  of  the  large  blocks,  the  four  columns  represent  the  4 
local  time  zones  (with  indices  11,  12,  21,  and  22)  while  the  4 
rows  represent  the  4  Kp  zones  in  order  of  increasing  magnetic 
activity. 

(4)  The  contents  of  any  cell  up  to  cell  495  of  December  (Key  =  11  x 
864  +  495  =  9999)  can  be  actually  displayed  by  using  the  EDIT 
processor  of  the  Honeywell  Sigma  9  computer.  For  greater  keys, 
one  must  resort  to  the  KFM  (keyed  file  manipulation)  processor. 

Having  accessed  at  most  48  cells  (4  Kp  x  4  local  time  x  3  months)  one 
must  first  adjust  the  data  to  the  user’s  sunspot  number,  magnetic  activity, 
and  month.  The  values  of  the  characteristics  will  actually  be  computed  for 
all  local  times  (hourly  intervals) . 

2.3.2  Correction  to  the  User's  Rz  and  Kp 

For  each  of  the  (at  most)  48  cells  the  foF2  is  adjusted  to  the  user's 
sunspot  number  Rz  using: 


foF2  -  a(Lt,  M,  Kp)  Rz  +  b(Lt,  M,  Kp) 


(2.10) 
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For  the  remaining  characteristics  (foF2  percentile,  median  spread, 
spread  percentile)  the  following  is  used: 

Uc  =  Sc  +  a  (Lt,  M,  Kp)  (Rz  -  Rz)  (2.11) 

where  Sc  is  the  characteristic  computed  from  the  standardized  histogram, 

Rz  is  the  average  sunspot  number  of  the  cell, 

Uc  is  the  characteristic  corrected  to  the  user's  Rz. 

As  discussed  above  a(Lt,  M,  Kp)  is  taken  to  be  zero  for  the  spread  (see 
Figure  8(a)). 


At  this  point,  for  each  of  the  12  (Lt,M)  pairs,  there  exists  an  array 
of  at  most  4  pairs  (Kp^,  Uc^),  i  =*  1(1)4,  where  Uc^  is  the  value  of  the 
characteristic  adjusted  to  the  user's  sunspot  number  (if  necessary).  To 
interpolate  to  the  user's  Kp,  UKp,  the  following  is  used: 


Ucl 

UKp  £  Kp 

P3(UC;l,  Uc2,  Uc3,  Uc4) 

Kp3  <  UKp  <  Kp 

Uc4 

UKp  >  Kp 

(2.12) 


where  is  the  Lagrange  interpolation  polynomial  of  degree  3  through  the  4 

points.  If  there  are  less  than  4  such  pairs  (due  to  non-existent  cells  or 
"missing"  Rz  coefficients  a,b)  the  degree  of  the  Lagrange  interpolation  is 
appropriately  reduced.  If  the  number  of  points  available  is  zero,  there  will 
be  less  than  12  triples  to  be  used  in  interpolation  in  the  next  section.  In 
such  a  case,  the  user  is  informed  and  the  interpolation  procedure  in  section 
2.3.3  is  omitted.  However,  the  user  can  obtain  the  available  triples  (<  12). 

2.3.3  Interpolation  to  Month  and  Local  Time 


The  12  triples  (M^,  Lt^,  Uc^),  i  =  1(1)12,  corresponding  to  the  6  half 
months  of  data  with  2  periods  of  local  time  each  are  used  in  a  least  squares 
model  to  obtain  an  array  of  the  characteristic  at  the  middle  of  the  user's 


month  and  for  each  hour  of  local  time, 
monthly  variation  is 

/  i 


Thus  it  is  essentially  assumed  the  data 
middle  of  the  half  month.  The  form  of 


The  variable  used  to  describe  the 
i  is  odd 
i  is  even 

for  a  given  half  month  occur  at  the 
the  model  to  fit  the  data  is: 


4 

Uc(Lt,W)  =  bn  +  Z  b.,_ . 

°  j-1  2j  1 


(simr  Lt/12)^  +  b2j  (cosir  Lt/12)^ 


+  l  b  (W/12)J 

j-1  j+8 


W  in  [1,  12]  ;  Lt  in  [0,  24] 


(2.13) 
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The  subroutines  STEP1  and  STEP2  perform  a  stepwise  regression  procedure 
for  calculating  and  selecting  the  statistically  significant  (non-zero) 
coefficients  bj:j=0(l)ll.  See  also  Ralston  and  Wilks  [24]  or  Frazer  [8], 
These  subroutines  require  a  pair  of  parameters  (F1,F2)  which  are  used  for 
entering  and  removing  a  given  variable  in  the  model  (2.13).  Typical  values 
of  FI  and  F2  are  about  5  for  a  95%  level  of  significance. 

2.3.4  Illustration 

A  display  of  the  results  of  the  main  prediction  program,  PREDSYS,  is 
now  given  where  the  user  has  input  the  conditions: 

Long  =  67° W,  Lattg  =  45°N,  Rz  =  30,  Kp  =  1+,  M  =  2 

He  also  desires  a  90%  -  tile  for  each  of  the  characteristics.  See  COMPUTER 
OUTPUT  2  for  an  example  of  the  actual  execution  procedure. 

In  COMPUTER  OUTPUT  3,  the  foF2  (MHz),  spread,  and  slab  thickness  (km) 
are  listed  for  each  hour  of  local  time.  Finally,  in  COMPUTER  OUTPUT  4,  the 
actual  cells  of  ALLMON  (directly  accessed  in  the  above  predictions)  are 
displayed.  See  remark  4  in  section  2.3.1. 

In  Table  7,  the  actual  (significant)  regression  coefficients  of  (2.13) 
as  well  as  the  Root  Mean  Square  Deviation  of  the  estimates  are  displayed, 
for  each  of  the  4  characteristics,  corresponding  to  the  set  of  predictions 
in  COMPUTER  OUTPUT  3.  The  empirical  model  for  slab  thickness  is  described 
in  Chapter  3. 

In  Figures  9(a),  9(b),  and  9(c)  contour  maps  of  foF2,  spread  and  slab 
thickness  are  displayed  for  the  above  input  conditions,  including  all 
northern  latitudes. 

In  conclusion,  the  following  remarks  are  made: 

(1)  The  Root  Mean  Square  Deviation  (in  Table  7)  is  only  strictly 
applicable  to  the  deviations  between  the  actual  (12)  observed 
values  of  a  characteristic  and  the  calculated  values  using  (2.13). 
However  these  deviations  should  approximately  apply  for  any 
predictions  Uc  (Lt,W)  in  view  of  the  natural  continuity  of  the 
characteristic,  in  local  time  and  month. 

(2)  The  accuracy  of  estimates  (and  even  the  capability  of  prediction) 
is  dependent  on  the  number  of  ionograms  available  in  the  various 
cells  directly  accessed.  Inspection  of  the  data  totals,  for  a 
typical  month,  (see  Table  4)  suggests  the  following  real  spatial 
domain,  D,  for  the  model. 

D  =  (A1  x  A2)  U  (B1  x  B2)  U  (Cl  x  C2)  (2.14) 

where  Al,  A2,  Bl,  B2,  Cl  and  C2  are  the  sets: 

A1  :  [ 180° W,  135°W] 

A2  :  [  45°N,  90°N] 


C 


A 


I 


Bl  :  [135°W,  90°W] 

B2  :  [  30°N,  90°N] 

Cl  :  [  90°W,  45°W] 

C2  :  [  0°  ,  90°N] 

The  first  set  of  each  pair,  in  (2.14),  is  the  longitude  range;  the 
second  is  the  geomagnetic  latitude  range. 

In  addition  the  sunspot  number  range  is  restricted  to  Rz  112 
(Table  3)  although  there  is  some  evidence  [6]  that  foF2  is 
independent  of  Rz  for  Rz  >  150.  A  more  serious  limitation  is  that 
approximately  two  thirds  of  the  data  is  in  the  low  sunspot  number 
range  (Table  2) 

(3)  One  could  directly  access  5  months  of  data  with  little  extra 
effort  by  redefining  the  size  of  the  arrays  in  PREDSYS.  This 
would  yield  at  most  20  triples  of  observations  to  be  used  in  the 


TABLE  7 

Illustration  of  Significant  Regression  Coefficients  in  Equation  (2. 13) 

COEFFICIENT  b. 
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COMPUTER  OUTPUT  Z  Example  of  Control  Procedure  Required  for  Predictions 


.EXAMPLE  SET  OF  PREDICTIONS  USING  EXPREDSYS 
!  .FOR  LONGITUDE =67  W(DEGREES) ,LATITUDE=45  N (DEGREES) 


!.RZ=30, FEBRUARY, KP=1+ 

!X£Q  EXPREDSYS ,36 
0036  -  i . 

0037  -  ! .A  FILE  NULL  IS  BUILT  CONSISTING  OF  A  SINGLE  RECORD. 
0038  -  !. NAMELY  THE  2  CHARACTERS  t. 

0039  -  ! .THIS  IS  USED  TO  BREAK  THE  EXECUTION  AFTER  THE 
0040  -  (.USER  OBTAINS  HIS  CONTOUR  DATA, ON  UNIT  80, TO 
0041  -  1 .BE  USED  FOR  EDITING  .  AFTER  HE  COMPLETES  THE 
0042  -  f. TASK, HE  CAN  RESUME  EXECUTION  AT  THE  LINE  91 
0043  -  ! .FOR  CONTOUR  PRODUCTION. 

0044  -  !B  NULL 
1.000  ! . 

2.000 
0045  -  ! . 

0046  -  i  .THE  OLD  GRID  OF  CONTOUR  DATA  IS  DELETED. 

0047  -  I . 

0048  -  ! DELETE  NEWPL0T1 
$ 

0300  FILE  DOES  NOT  EXIST. 

0049  -  1SET  F:105  ME;IN 
0050  -  (SET  F : 6/ALLMON ; IN ;KEYED;DIRECT 
0051  -  1SET  F:7  ME;OUT 
0052  -  ISET  F : 80/NEWPLOT 1 ;OUT ;KEYED;DIRECT 
0053  -  (SET  F:70/ADJRZ;0UT 
0054  -  !SET  F:71/ADJKP;0UT 
0055  -  ISET  F : 108/REGFIT;0UT 
0056  -  ILYNX  BPRED6YS , BMISC , BSSTEP 1 , BSSTEP2 
:P1  ASSOCIATED. 

GR0UPC2  OF  : LIB. -SYS  ASSOCIATED  FOR  POSREC 
GROUP 04  OF  : LIB.: SYS  ASSOCIATED  FOR  SQRTF 
GR0UP21  OF  :LIB. :SYS  ASSOCIATED  FOR  9DEFILE 
GR0UP23  OF  : LIB.: SYS  ASSOCIATED  FOR  7ERRARG 
GROUP  14  OF  :LIB.:SYS  ASSOCIATED  FOR  9IFR 
GROUP  17  OF  :LIB. :SYS  ASSOCIATED  FOR  9IMPUTL 
•  •  ALLOCATION  SUMMARY  •  • 

PROTECTION  LOCATION  PAGES 

DATA  (00)  A 000  9 

PROCEDURE  (01)  B600  10 

DCB  (10)  B200  2 

0057  -  I. 

0058  -  !  .ST  GIVES  CURRENT  CPU  TIME  CONSUMED. 

0059  -  I. 


COMPUTER  OUTPUT  3.  Example  of  Predictions  Using  Control  Procedure  of  COMPUTER  OUTPUT  2 


CPU  :  .  CON  =  00:  19:0-4  [NT  =  25  CHG  :  $  .94  CALS  =  1255 
0061  -  !S 

SPECIFY  YOUR  LOCATION _ EXAMPLE, FOR  OTTAWA - 

-75. 0,M5.,  56, 2<  1.7-,.  95,.  75,...  MEANS 
GEOGRAPHIC  LONGITUDE  IS  7SW ( DEGREES ) . 

GEOGRAPHIC  LATITUDE  IS  45N (DEGREES) . 

SUNSPOT  NUMBER  IS  56... IF  SUNSPOT  NUMBER  IS 
>  150, ENTER  150  AS  THE  CHARACTERISTICS  ARE 
ASSUMED  INDEPENDENT  OF  IT  IN  THAT  DOMAIN. 

MONTH  IS  FEBRUARY, IE  MIDPOINT  OF  THE  MONTH 
3-HOUR LY  KP  IS  2-. 

A  95  %-TILE  OF  F0F2  IS  DESIRED. 

A  75  %-TILE  OF  SPREAD  IS  DESIRED. 

ALL  CHARACTERISTICS  ARE  COMPUTED  FOR  EACH 
HOUR  OF  LOCAL  TIME. .  .INCIDENTALLY, A  95%-TILE 
IS  THE  VALUE  OF  THE  CHARACTERISTIC  SUCH  THAT  . 

95  %  OF  THE  DATA  IS  LESS  THAN  THAT  VALUE. 

FOFP  IS  IN  MHZ... SLAB  THICK  IS  IN  KM 
TO  STOP, ENTER  1000,0,0,0,0,0,0. 

ALSO  USE  THIS  AS  THE  LAST  RECORD  OF  AN  INPUT  FILE 
OF  CONDITIONS... IE  WHEN  PREPARING  CONTOUR  MAPS. 

SEE  COMMENT  STATEMENTS  IN  SOURCE  PROGRAM  ABOUT 
PREPARING  CONTOUR  MAPS. 

?-67,45, 30,2, 1.3, .9, .9  . 

THERE  ARE  12  TRIPLES  OF  IDCAL  TIME  AND  MONTH  DATA 
IF  LESS  THAN  12, TRY  AN  ADJACENT  CELL. .  .OR  DIRECTLY 
ACCESS  THE  CELLS  USING  EDIT  OR  KFM  PROCESSORS. 

THE  PRELIMINARY  INTERPOLANTS  ARE  ON  UNITS  70  &  71 
GEOGLONG:  -67.0  GEOGLAT:  45.0  SUNSPOT:  30. 

MONTH:  2  KP=1.3 

PROB  F0F2:  .90  PROB  SPREAD:  .90 

GEOMATLAT:56.4 


LOCAL  TIME 

F0F2 

F0F2% 

SPREAD 

SPREAD% 

SLAB  THICK 

0 

1.5 

3.0 

2.0 

2.5 

181. 

1 

1.6 

3.0 

2.0 

2.5 

175. 

2 

1.6 

3.0 

2.0 

2.6 

178. 

3 

1.6 

2.9 

2. 1 

2.7 

193. 

4 

1.6 

2.8 

2.1 

2.8 

213- 

5 

1.7 

2.8 

2.1 

2.9 

231. 

6 

2. 1 

3.1 

2.1 

2.9 

238. 

7 

2.8 

3. '8 

2.1 

2.9 

231. 

8 

3.7 

4.7 

2.1 

2.8 

213- 

9 

4.6 

5.6 

2.0 

2.7 

193. 

10 

5.3 

6.3 

1.9 

2.6 

178. 

11 

5.6 

6.7 

1.9 

2.5 

175. 

12 

5.7  ‘ 

6.8 

1.8 

2.5 

181. 

13 

5.7 

6.7 

1.9 

2.5 

192. 

14 

5.7 

6.7 

1.9 

2.6 

201. 

15 

5.7 

6.8 

2.0 

2.7 

202. 

16 

5.7 

7.0 

2.1 

2.8 

198. 

17 

5.6 

6.9 

2. 1 

2.9 

192. 

18 

5.2 

6.6 

2. 1 

2.9 

189. 

19 

4.5 

5.9 

2.1 

2.9 

192. 

20 

3.6 

5.0 

2.1 

2.8 

198. 

21 

2,7 

4.1 

2.1 

2.7 

202. 

22 

2.0 

3.4 

2.0 

2.6 

201. 

23 

1.6 

3.1 

2.0 

2.5 

192. 

24 

1.5 

3.0 

2.0 

2.5 

181. 

71000,0,0,0,0,0,0 
•STOP*  0 


COMPUTER  OUTPUT  4.  Extract  From  the  File  ALLMON 


fi.fi 

£i 


3  ' 
& 


J 

*-> 

S 

! 


8 

sl 


< 

ui 

a 

u 


8R. 


nnHRMN^oinnNNNM 

eoe-o»-Nwooineo 

mNnoNNonn»>«r<<oN 
•  ••••••••••»»• 

GD")«<-eDai«-n»>n»>\BtD 

xN4nNNOIPO<04feV<* 

MnNomiotnoncDKinoir) 
•  «•••••••••••• 

irj’onrv^^Morvio^^rirt^N 

<N<mop<no-o»fVNf) 

^noKNO'Orses^ooiO(^». 

NN««ODaCOOO<N4 

N  NNN  CM  CM  CM 

WNNONMOOW-ONMir 
inf  noMnNN'CnniM'io 

CM  ▼  FV  CM  ▼  ID  CM  *  IN  V 

<«M-nnNirinNnon< 
«  yvm  oeaaiNMM  o»>  » 

▼  aina<^in(Nryrs^^r>o 

N  NNN  N  NN  N  N  N  N  Nlv  N 

imoriimniDnwinin  rinnci 
CM^n^einNin^MCMtNrMCM 

«N»>O'®«'0ONKWnifl 
NNNMNNNNNKMNNN 
1  I  It  1111-11  till 


I  mom  m 

*  ««  r>  n  •< 


n  cv 

«  m 


NNNNI 

«  w  n  «  f 


tnt  **cm  n«»Mn 

INN-r-nNN 
m  «*  «*  CM  CM  CM  IM  CM  CM 

^NNNNNKKN 

tpmnnnnnm 


-  Z  M>  C  O  ©  © 
«  £  i  O  ©  O  © 
£•>0000 
►  J  N  •  •  •  • 

-w<« MUM 
3  C  ►  N  N  N  N 

U.  I  »-•«-«•«  *o 


§§< 
o  o  < 


It 


^8 

!«<« 


58 


n  « 
a  cd 
+  © 


«•  ft 
CD  CD 
<•  « 


£  £ 


IN«N 

CM 

I  <►  ft  ft 


n  cm  « 
o  «  < 


•“-M 

CD©© 


«ft 


»  Bfr> 

w  w® 


i*«oo«oo 

I  vt  M« 

INmmNmm 

N-O-NtN 

©  CD  ©  w  +  ▼ 

o  ®  ►  r>  ©  «r 
DOMNMf 


.  ft  ©  m  r 
i  «  W  d*  < 


in  n  io«  — 

«.ft  O  OB  CM  O 
1  :«  »  ftft  W  © 


iinwownwmoipcm© 
(CjMMBKMnN®  win 


n  '0MM»*8fM,n-««npiMnN 
•  •*«••••••••••• 

ft  ftftftftfvBftftftftftftftft 

ncnommncnnnnr«in«c» 
ft  »nm^d>»-db«n«m»o 


i  ®  —  cm 

-CNN 
I  11 


MOTN 
<  © ©• 


m  •>  n  ©  *  cm  o  ®  e  n  o 

KNN-C-C-OMnNNN 
11  1  11  1*1*1  1  1  t 


I  ftCMCMftCMCMftftftftftCMCMW 

*  npim -nn  M*<NnM-Nf> 
i  ***««NNMnf>«<M«rtNN 
■  ftMMMMMMNNNNNflN 
.  NNNKNNNNNNNKNN 

i  nnBnnnnnfinrtinnn 


piMnooooNNf  om  n 

w» 

N-.wfNn-rB«MwaN 

M  ft  ©  ft  ©  ft  N  *0B  Cl  •"  M  ©ft 
~  ©  »  Iff  **  fttt©©CNW©»-ft 
•  ••••••••••••• 

tMKN^IHf  rip-OMiunMin 


n<r o««no< 
tiff  in  *  *«AMMIDCV 


ft©ft©nmcM©ftft©©©n 

®©©o©in©*M~'©inrtfv 

fM  *  ©  CM  n  ft©©  CM  * 


ft©M-®©©ft©ft©*®ft© 


okvbdn  nD>8«<ctrnn 

a*  a  ft  »  ft  in  mnnbdm 

■CMN  CNN  ft©©©ftftftft 
I  t  t  I  I  t  1  I  I  t  1  t  1  1 


©  © w  -«  n  ©  T<e«innBn 

n«Bne«  non(M»D>».n 


nnnnmnnn  ft  nnnnn 
-Nwc  wnn-Miiwrin 

«  mm  m  m  n  rtMMMMMtn 
MMMMMMMnrMCMCKMtin 
N  NN  N  NN  rvryNrs.NN.NN 

tnrtrjnnn  r»n««r>wnn 


iiiiiii  is  iiiiiiiiiiiiii!  a  !i§if§i§iii§§; 


Oft© 

i  rt  w© 

i  winin 


«•  ft  W  t«<NB^O- 


—  cm  r>  ©  ri  ©  ft*o«Nn©in 

©©©©©©  OOmmmmmm 
Iff tit fMMMMMMM 

CM  CM  CM  CM  CM  CM  CMftftCMftCMftftD 


6  12  18 
GEOGRAPHIC  LOCAL  TIME  (hrs.) 


Figure  9(a).  Contours  of  Characteristics  in  Latitude  and  Local  Time,  foF2  (MHz). 
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6  12 
GEOGRAPHIC  LOCAL  TIME  (hr*.) 

Figure  9(b).  Contours  of  Characteristics  in  Latitude  and  Local  Time,  spread  in  f  1,4/. 
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Figure  9(c).  Contours  of  Characteristics  in  Latitude  and  Local  Time,  Slab  Thickness  (km). 


determination  of  the  model  (2.13).  Balanced  against  this,  there 
would  be  more  monthly  variation  to  be  described  by  the  model. 

(4)  The  current  cell  alignment,  for  a  given  month,  leads  to  a  prolific 
number  of  cells  with  few  ionograms  (especially  for  those  corres¬ 
ponding  to  high  magnetic  activity).  Perhaps  only  2  Kp  zones,  or 
a  single  zone  for  all  Kp's  would  ameliorate  this  situation.  A 
"smoothed  bi-weekly  Kp  average"  could  be  attached  to  each  of  the 
resulting  216  (=  864/4)  cells.  To  consider  such  questions  one 
only  needs  to  apply  ACCMONB96  which  operates  on  the  file  MONB96 
(blocked  cells,  but  unstandardized) .  These  files  have  been  stored 
on  magnetic  tape  (all  months);  see  Appendix  A  (listing  of  tapes). 


3.  N(h)  MODEL 

3.1  DESCRIPTION  OF  THE  DATA  BASE 

About  80,000  Alouette  1  and  2  ionograms  have  been  scaled  to  produce 
almost  complete  topside  electron  density  profiles.  These  profiles  have  been 


35 


stored  on  16  9-track  tapes.  Alouette  2  was  launched  in  November  1965,  and 
unlike  Alouette  1  it  had  an  elliptical  orbit:  apogee  =  3000  km,  perigee  = 

500  km.  Since  the  ionograms  were  scaled  on  demand  they  were  actually  stored 
on  tape  as  the  need  arose,  resulting  in  a  very  fragmented  data  base,  unlike 
ALOSYN.  Indeed,  the  actual  details  of  the  reduction  procedure  [11],  to 
generate  N(h)  profiles,  was  varied  throughout  the  years. 

In  Table  8  a  summary  of  the  data  totals  of  the  various  tapes  over  the 
years  is  presented.  Some  remarks  pertaining  to  this  table  follow. 

(1)  The  ionograms  within  a  volume  are  in  chronological  order  (as  they 
were  scaled) . 

(2)  In  A  Format,  the  electron  density  profiles  have  been  calculated 
by  assuming  fN,  the  plasma-frequency,  is  linear  in  height  on  each 
lamination. 

(3)  In  B  Format,  the  reduction  has  been  achieved  by  assuming  the  log 
of  the  electron  density  varies  parabolically  in  each  lamination 
(but  linearly  in  the  first  lamination  nearest  the  satellite). 

(4)  The  profiles  tabulated  by  the  Ames  Research  Centres  (AMES)  are 
actually  interpolated  to  give  electron  densities  at  every  one 
hundred  kilometres  from  the  height  of  the  satellite.  In  addition 
they  distinguish  between  field-aligned  and  non  field-aligned 
profiles. 

(5)  The  topside  TEC,  in  A  Format,  was  calculated  by  integrating 
(exactly)  the  model  electron  density,  N(h),  from  the  height  of 
the  satellite  to  a  height  of  maximum  electron  density.  The  latter 
was  calculated  by  passing  a  parabola  through  the  lowest  2  scaled 
points  and  assuming  foF2  was  0.01  MHz  greater  than  the  lowest 
scaled  fN.  Thus,  the  topside  TEC  has  been  adjusted  to  include 
only  the  area  to  the  last  scaled  point;  the  remainder  of  the  TEC 
(bottomside)  will  be  computed  from  the  model  to  be  described  in 
section  3.2.2.  In  B  Format,  the  topside  TEC  was  computed  to  the 
lowest  scaled  point,  assuming  the  electron  density  is  linear  on 
each  lamination.  The  Ames  profiles  do  not  give  TEC  so  a  value 
consistent  with  the  B  Format  profiles  was  calculated  and  put  on 
tape. 

There  also  exist  10  tape  volumes  of  ISIS  1  and  2  profiles  in  B  Format. 
These  two  satellites  are  similar  to  Alouette  2  and  1,  respectively,  in 
orbital  configuration.  They  were  launched  in  1969  and  1971  during  a 
decreasing  period  of  solar  activity.  Thus,  the  data  reduction  procedures  to 
be  described  below  are  readily  adaptable  to  these  data  as  well. 


3.2  DEVELOPMENT  OF  THE  MODEL 

3.2.1  Organization  of  the  Data  Base 

The  intention  of  this  section  is  the  construction  of  a  chronological 
ordered  data  base  of  the  N(h)  profiles.  Each  profile  of  this  new  base  will 


36 


TABLE  8 

Nth)  Profile  Data  Available 


TAPE 

VOLUME 

SATELLITE 

FORMAT 

ERA 

NUMBER  OF 
I0N0GRAMS 

1  -  4 

ALOUETTE  1 

A 

1962  -  66 

43700 

5  -  7 

ALOUETTE  1 

B 

1963  -  70 

18400 

1  -  5 

ALOUETTE  2 

B 

1965  -  72 

12000 

6-9 

ALOUETTE  2 

AMES 

1966  -  72 

5200 

1  -  5 

ISIS  1 

B 

1969  -  present 

35000  ‘ 

1  -  5 

ISIS  2 

B 

1971  -  78 

43000 

consist  of  the  identification  parameters  (time,  location,  Kp,  and  Rz) 
together  with  certain  profile  characteristics,  e.g.  TEC  (topside),  TEC 
(bottomside) ,  hmax,  and  slab  thickness.  The  definition  of  these  characteris¬ 
tics  will  be  considered  in  more  detail  in  the  next  section.  This  new  data 
base,  ALOUETTE  NH,  will  then  be  similar  to  ALOSYN,  although  not  nearly  as 
extensive  of  course.  Thus,  the  entire  procedure  of  Chapter  2,  with  a 
necessarily  coarser  grid,  can  be  applied  to  ALOUETTE  NH. 

An  illustration  of  this  approach  will  be  given  in  section  3.2.4,  for  a 
selection  of  10,000  Alouette  1  profiles  (1965-70)  which  cover  (sparsely)  a 
half  solar  cycle,  i.e.  Rz  in  [12,  112]  (see  Table  3).  The  empirical  model 
used  for  slab  thickness,  in  COMPUTER  OUTPUT  2,  will  also  be  fitted  to  these 
data. 


The  process  of  constructing  ALOUETTE  NH  consisted  of  processing 
(calculating  the  frequency  and  heights,  calculating  the  topside  TEC,  and 
attaching  the  Kp)  for  each  of  the  16  volumes,  described  above,  one  at  a  time. 
The  resulting  16  output  files  were  concatenated,  sorted  chronologically,  and 
blocked  (at  the  rate  of  10  profiles  per  record).  Only  profiles  in  the 
longitude  range  [-180°,  -45° ]  are  included  in  ALOUETTE  NH.  However,  due  to 
the  small  data  base,  all  geomagnetic  latitudes  are  included,  and  it  is 
assumed  that  slab  thickness  is  symmetric  about  the  geomagnetic  equator. 

The  data  actually  kept  for  a  particular  profile  of  ALOUETTE  NH  are  as 
follows : 

Volume  (tape  number  of  original  profile  of  ALOUETTE  N(h);  see 
Table  8) 


Kp  (integer  part  only) 

Dip  (dip  angle  at  the  satellite,  degress) 

Job  (see  below) 

Year,  Day,  Hour,  Minute,  Second  (Universal  time) 

TTEC  (topside  TEC  from  hQ  to  in  electrons /cm2,  where  hQ  is  the 
satellite  height  and  hL  is  the  lowest  height  of  the  topside 
profile) 

L+l  (number  of  pairs  of  points  for  profile) 

fH  (gyrofrequency  at  hQ,  MHz) 

Long  (degrees) 

Lattg  (degrees) 

X  (Solar  angle  at  satellite,  degrees) 

fN  (a  4x1  vector,  plasma  frequency  at  h 

MHz)  ° 

h  (a  4x1  vector,  real  height  hQ,  ^  2*  ^  1*  and  hL*  km) 

The  Job  parameter  is  a  subjective  index  concerning  the  quality  of  the 
profile  on  a  scale  of  4  (excellent)  to  9  (scaled  under  protest).  In  the 
case  of  the  Ames  profiles,  this  parameter  was  used  here  to  indicate  whether 
the  profile  is  field  aligned  or  not, 

0  ,  non  field-aligned 

Job  = 

Year,  field-aligned 

In  the  field-aligned  case,  the  longitude  and  latitude  are  the  arithmetic 
means  of  the  longitudes  and  latitudes  of  the  various  points  along  the  magnetic 
field  direction. 

3.2.2  Definition  of  the  Profile  Using  ALOUETTE  NH 

Here,  a  method  of  defining  a  parabola  will  be  outlined,  using  the  data 
of  ALOUETTE  NH,  to  represent  the  profile  from  hL  down  to  h.  ,  a  base  for  a 
hypothetical  F2  layer  to  be  defined  below.  From  this  parabola,  rough 
estimates  of  the  bottomside  total  electron  content  (BTEC),  hmax,  and  slab 
thickness  can  be  made.  High  quality  estimates  will  also  be  compared  with  an 
alternative  method  (for  limited  samples  of  data)  using  ground  echoes.  High 
quality  ground  echoes  usually  only  occur  on  the  Alouette  ionograms  from 
soundings  taken  over  water  or  over  less  densely  populated  areas  where  noise 
levels  are  low.  This  fact,  combined  with  the  additional  scaling  effort  that 
would  be  required,  precludes  the  use  of  this  method  for  the  entire  data  base. 


»  **L-2»  hL-l*  and  V 
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The  first  step  is  to  define  certain  subsets  of  ALOUETTE  NH  consisting 
of  "quality"  profiles.  Here,  quality  means  that  the  last  point  (NL,hL)  is 
"reasonably  close"  to  the  peak  electron  density  (Nmax,  hmax) .  The  electron 
density  N  is' related  to  the  plasma  frequency  f^  at  any  height,  h,  by 

N  =  1.24  x  104  x  f ?.  (3.1) 

N 

where  N  is  in  electrons/cm3,  f^  is  in  MHz.  Usually  the  profiles  from  topside 
soundings  tend  to  be  too  low  [12],  However,  Jackson  [12]  notes  that  they  can 
also  be  high  (particularly  for  soundings  when  the  satellite  is  above  2000  km) . 
To  define  the  subset  of  quality  profiles  some  weak  conditions  will  be  imposed 
on  the  pairs  of  points,  rather  than  simply  relying  on  the  subjective  job  index. 
Before  enumerating  these  restrictions,  a  comparison  of  the  fN  at  h=hL  and 
foF2,  obtained  from  ALOSYN,  for  some  selected  ionograms  is  given. 

In  Table  9,  A  is  defined  as: 

A  =  100  [  foF2  -  W]/W  (3‘2) 

From  this  Table,  the  median  A  is  between  0  and  2,  i.e.  foF2  exceeds 
f„(hL)  by  2%  on  average.  Also  about  85%  of  the  data  are  such  that  A  <  5. 

These  profiles  should  comprise  the  quality  subset. 

TABLE  9 

Comparison  of  Lowest  Scaled  N(h)  Profile  Point  and  foF2  From  ALOSYN  for  Selected  Passes 


A 

CUMULATIVE  PROBABILITY 

DISTRIBUTION  OF  A 

FEBRUARY 

APRIL 

FEBRUARY 

APRIL 

1965 

1965 

1967 

1967 

R  z  --  1? 

Rz=  1? 

Rz  =  76 

Rz=  81 

-6 

4 

1 

7 

1 

-5 

5 

? 

8 

3 

-4 

6 

3 

9 

3 

-3 

7 

5 

10 

3 

■? 

11 

6 

16 

3 

-i 

IB 

11 

21 

11 

0 

?5 

16 

26 

27 

1 

51 

37 

53 

61 

? 

65 

5? 

61 

83 

3 

/6 

68 

73 

89 

4 

8? 

79 

78 

92 

5 

88 

88 

83 

96 

6 

91 

90 

86 

96 

i  ’ 

93 

9? 

87 

97 

n 

94 

93 

89 

97 

0 

95 

96 

90 

100 

in 

96 

96 

90 

100 

NIIM8IR  Of 

IONOGRAMS 

589 

275 

763 

71 

.  . . 
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The  weak  restrictions  for  defining  the  quality  subset  of  ALOUETTE  NH 
are  now  enumerated: 

ho  N  +Nt 

<»  WV  i  /  N(h)dh  i  -V1  <W 

(2>  nl-2  '  "l-1  '  \  >  \-2  ’  hL-l  *  h, 

(3)  h  <  2000 

o 

(4)  <  600 

Foi.  the  base  of  the  layer,  h^,  at  which  N^O,  the  following  heuristic 
model  based  on  experience  [E.L.  Hagg,  private  communications,  1978]  has  been 
used: 

2.5  Lt  +  230  ,  Lt  in  [0,4] 

30sin(7rLt/8)  +  210  ,  Lt  in  [4,12]  (3.3) 

180  exp  [0.0204(Lt-12) ]  ,  Lt  in  [12,24] 

where  Lt  is  local  time  in  hours,  and  h^  is  in  km. 

An  empirical  model  for  the  height  of  maximum  electron  density  of  the  FI 
layer  [25],  given  by 

hmax  FI  =  165  +  0.6428  x»  (3.4) 

is  compared  with  (3.3)  in  Figure  10  (x  is  in  degrees). 

The  complete  density  profile  is  summarized  in  Figure  11.  A  few  remarks 
are  added: 

(1)  The  following  artifice  will  be  convenient.  For  the  exponential 
function 

N  -  Nx  expKh-h^/H]  (3.5) 

passing  through  any  2  points  (N2»h2)  and  (Ni.hj)  (with  N  in  [N2.N1] 
and  h  in  [h  ,h2]), 

H  “  (hl  '  h2>/iln(W  <  0  (3.6) 

is  called  effective  scale  height.  Note  that  the  symbols  in  (3.5) 
and  (3.6)  are  generic;  they  will  be  specified  in  (2)  and  (3) 
below. 


(2)  In  a  few  cases,  i.e.  for  Alouette  2,  the  satellite  is  considerably 
below  1000  km.  In  such  cases,  an  exponential  curve  is  passed 


HEIGHT  (km) 
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- HEIGHT  OF  MAXIMUM  DENSITY  OF 

Fi  LAYER  (hmax.  FI) 

JUNE,  GEOGRAPHIC  LATITUDE  40°N 


HEIGHT 


i - 1 - 1 - r 

No  Nl_2  Nl_,  Nl 
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through  the  points  (N0,h0)  and  N^_2»  h^_2)>  in  order  to  make  a 
correction  for  the  total  electron  content  above  the  satellite. 
This  extra  area  is: 


ETEC 


QO 

-I 


N(h)dh  -  -N  H 
o 


(3.7) 


(3)  The  constants  of  the  parabola  are  chosen  in  order  that  it  is 

continuous  and  has  a  continuous  slope  at  (Nn^).  The  slope  of 

N  -  Nl  exp  [(h-h^/Hj 

at  h  *  hL  is  used  for  this  purpose,  where  H  is  a  scale  height  for 
the  exponential  on  [hL,  hL_]_].  T^e  requirement  is  that 

N(hb)  =  0  where  hjj  is  the  base  of  the  layer,  as  previously  defined. 
Thus  the  parabola,  illustrated  in  Figure  11,  is 


N  =  N_ 


Setting  N'(h)  =  0,  in  (3.8) 

hmax  ■  h^  -  (b/2c)  (note  b,c  <  0  since  H  <  0  and  h^  <  h^) 
Nmax  *  a  +  b(hmax  -  ly)  +  c(hmax  -  1^) 


Finally 


BTEC 


N 


(h)dh 


(3.8) 

(3.9) 

(3.10) 

(3.11) 


and  slab  thickness  is 

Th  -  (ETEC  +  TTEC  +  BTEC) /Nmax 

3.2.3  Definition  of  Profile  Using  Ground  Echoes 


(3.12) 


In  this  section,  an  alternative  approach  for  estimating  profile 
characteristics  by  adapting  an  idea  of  Dyson  [7]  to  the  present  data  base 
will  now  be  outlined.  The  method  essentially  is  an  attempt  to  model  a 
complete  profile  (topside  and  bottomside)  using  only  data  available  from 
satellite  soundings. 


For  frequencies,  fx,  at  which  ground  echoes  occur  the  group  path  due 
to  ionization  below  h^  is  given  by  (see  Figure  12) 
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FREQUENCY,  fx 


Figure  12.  Ground  Echoes  Used  for  Estimation  of  Bottomside  Total  Electron  Content 


?L 


p'  -  p; 

m  t 


=  j  (p*  -  i)  dh  + 


(3.13) 


where  p'  is  the  scaled  group  path  (virtual  height  or  apparent  range); 
m 

p^  is  the  group  path  from  hQ  to  h^; 

p'  is  the  group  refractive  index;  a  function  of  fx,  gyrofrequency 

f  ,  dip  angle,  and  electron  density,  N. 

H 

The  integral  in  (3.13)  is  transformed  by 


X  =  Xmax  | 1  -  [ (h-hmax) / 2T ] 2 } 

where  X  =  (fM/fx)2  =  N/12400fx2 
N 

T  =  quarter  thickness  of  the  parabola  to  be  determined. 
This  gives  the  single  linear  equation 


(3.14) 


(3.15) 


where  R 
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The  procedure  adopted  here,  for  the  solution  of  (3.15)  for  T,  is  as 
follows: 

(1)  A  number,  Ng,  of  pairs  (p^jfx)  are  scaled  and  smoothed.  Only  those 
pairs  satisfying  both 

(a)  fxF2/fx  <  0.85 

(b)  p'  -  h  >60  (km) 

m  o 

are  actually  used  in  the  calculation. 

(2)  The  fxF2  (and  foF2)  are  obtained  from  the  ALOSYN  data  base. 

(3)  The  integral  p£  is  computed  using  the  original  complete  profiles 
of  ALOUETTE  N(h),  assuming  p*  is  linear  in  each  height  lamination. 

(4)  The  integrals  on  the  right  side  of  (3.15)  are  computed  using  the 
obvious  trigonometric  substitutions  (to  remove  the  singularity) 
followed  by  an  application  of  the  composite  (closed)  5  point 
Newton-Cotes  integration  rule  [24] . 

(5)  For  each  of  the  Ng  pairs  of  scaled  echoes  determined  in  (1)  above, 
there  are  a  like  number  of  linear  equations: 

Rt  =  TAt  i  -  1(1)  Ng  (3.16) 

where  and  A^  are  respective  values  of  R  and  the  integrals  in 
(3.15).  Equation  (3.16)  is  solved  for  T  by  the  method  of  least 
squares 

N 

8  2 

i.e.  minimize  E  (R  -  TA  ) 

i-1  1  1 

Thus  the  parabola  (3.14)  on  [h^,  h^]  is  completely  determined,  so 

that  the  method  of  section  3.2.2  may  be  used  to  determine  an 

alternative  version  of  slab  thickness  (hmax  is  found  by  making 

(3.14)  continuous  at  (N  ,  hT)). 

L  L 

3.2.4  Illustration 


In  Tables  10(a)-(d)  and  ll(a)-(d),  slab  thickness  has  been  displayed 
for  various  sets  of  input  conditions.  These  values  have  been  computed  by 
the  procedure  outlined  in  section  3.2.2  and  then  sorted  into  a  (crude) 
partition  of  the  domain  of  the  ALOSYN  model,  see  below.  The  standard  error, 
Se,  of  the  mean  slab  thickness  is  given  by 

Se  -  S/N*  (3.17) 

where  S  »  standard  deviation  of  slab  thickness. 


N  =  ionogram  count 


Also,  the  foF2  value  displayed  Is  a  "check  value"  based  on  a  peak  electron 
density,  Umax,  which  is  the  ratio  of  the  (average)  TEC  for  the  cell  and  the 
(average)  slab  thickness. 


The  data  used  in  the  preparation  of  these  tables  consisted  of  approxi¬ 
mately  10,000  Alouette  1  profiles  in  the  time  period  1965-1970,  or  essentially 
half  of  the  sunspot  number  cycle,  Rz  in  [12,  112], 

The  domain  was  partitioned  into  384  cells  which  consist  of: 

(a)  2  longitude  zones  of  70°  each,  starting  at  180°W 

(b)  3  geomagnetic  latitude  zones  of  30°  each,  starting  at  90°N 

(c)  2  magnetic  activity  zones,  divided  at  Kp  =  2.5  (between  2+  and  3-) 

(d)  2  sunspot  number  zones,  divided  at  Rz  =  40 

(e)  4  local  time  zones  of  6  hours  each,  starting  at  midnight 

(f)  4  season  zones  of  3  months  each,  starting  at  January. 

In  Tables  10  and  11,  only  the  results  from  the  two  most  northern  geo¬ 
magnetic  latitude  zones,  and  for  the  eastern  longitude  zone,  have  been 
displayed  —  these  actually  contained  approximately  75%  of  the  profiles  used. 


TABLE  10(a) 

Mean  and  standard  deviation  of  stab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC, 
ionogram  count  and  foF2  (MHz)  for  Lattm  >  6(f  N,  Kp  <  2.5,  Rz  <  40. 


— — - 

\  Local 
^S\T  i  mo 

Month^\^ 

'0.6| 

16.121 

112.181 

118,241 

January 

to 

March 

232  77 

.25  96 

2.  3 

20R  70 

.?:■>  10R 

3.3 

219  32 

.32  211 

5.1 

226  50 

.34  188 

2.7 

Apr  i  1 

to 

June 

215  36 

.11  262 

3.6 

203  30 

.1?  435 

4.3 

248  32 

.20  721 

4.9 

230  33 

.21  327 

4.0 

July 

to 

September 

* 

- 

- 

- 

0c  tobnr 

to 

December 

- 

146  23 

.09  4 

5.7 

- 

- 

*  No  N(h)  profiles  available  (in  this  sample). 


TABLE  10(b) 


Mean  and  standard  deviation  of  stab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC, 
ionogram  count  and  foF2  (MHz)  for  Lattm  >  6(f  N,  Kp  >  2.5,  Rz  <  40. 


'  Local 

^sNTime 

[0,61 

[6.121 

[12,18] 

[18,24] 

Hcnth^'v\^ 

January 

to 

March 


261  40 

.39  4 

3.7 


287  27 

.10  6 

3  1 


4 

53 

5 

45 

July 

to 

September 


October 

to 

December 


TABLE  10(c) 

Mean  and  standard  deviation  of  slab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC, 
ionogram  count  and  foF2  (MHz)  for  Lattm  >  6(f  N,  Kp  <  2.5,  Rz  >  40. 


n.  Local 
^\Tlme 

[0,6] 

[6,12] 

[12,18] 

[18,24] 

Month 

January 

to 

March 


237  30 

.48  58 

3.  3 


250  25  281  38  285 

.15  127  .16  234  .28 

4.0  4.2  ' 


July 

to 

September 


October 

to 

Oecember 


323 

42 

.28 

7.2 

14 

285 

27 

.28 

5.4 

245 

260 

18 

.27 

5.9 

22 

237 

57 

.51 

7.0 

21 

0 

26 

6 

22 

4.4 

TABLE  10(d) 

Mean  and  standard  deviation  of  slab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC , 
ionogram  count  and  foF2  (MHz)  for  Lattm  >  60°  N,  Kp  >  2.5,  Rz  >  40. 


Local 

^\Time 
Month  Nv 

CO.  6] 

[6.121 

[12,18] 

[18,24] 

January 

to 

March 


July 

to 

September 


October 

to 

December 


311  71 

.41  15 

3.1 


377  100  ?8 

.34  21  .? 

3.? 


39  ?91  27  305 

11  .17  51  .29 


.28  10  .59 

6.8 


390 

16 

.42 

5.1 

9 

314 

41 

.29 

4.6 

95 

305 

32 

.29 

4.9 

16 

220 

47 

.59 

7.9 

16 

TABLE  11(a) 

Mean  and  standard  deviation  of  slab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC 
ionogram  count  and  foF2  (MHz)  for  Lattm  in  [30°  N,  6(f  N\  Kp  <  2.5,  Rz  <  40. 


Local 

^sTIme 

Month 

[0,6] 

[5,12] 

[12.18] 

[18,24] 

January 

to 

March 


July 

to 

September 


October 

to 

Oec ember 


!  41 

)  161 

3.5 

187 

.30 

A 

44 

452 

.8 

195 

.31 

5, 

26 

204 

,7 

,  57 

179 

31 

226 

35 

322 

.21 

1R8 

.53 

590 

2.7 

4 

.8 

6 

.0 

141 

4 

.20 

5.6 

4 

TABLE  11(b) 

Mean  and  standard  deviation  of  slab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC, 
ionogram  count  and  foF2  (MHz)  for  Lattm  in  [ 3(f  N,  6(f  N\,  Kp  >  2.5,  Rz  <  40. 


\.  Local 

^^Tlme 
Montn  N. 

[0.6] 

[6.12] 

'  [12.18] 

[18.24] 

January 

to 

March 


199  47 

.35  160 

5.6 


262  63 

.66  55 

2.9 


259  51 

.41  23 

2.7 


261  36 

.28  29 

4.4 


243  2 

.61  4 

5.4 


222  26 

.73  61 

4.0 


July 

to 

September 


October 

to 

December 


TABLE  11(c) 

Mean  and  standard  deviation  of  slab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC. 
ionogram  count  and  foF2  (MHz)  for  Lattm  in  \3(f  N,  6(f  N\  Kp  <  2.5,  Rz  >  40. 


Local 

^'v^Tlme 
Month  N. 

[0.61 

[6,12] 

[12.18] 

[18,241 

January 

to 

March 


July 

to 

September 


October 

to 

December 


410  13 

1.2  3 

1.4 


217  6 

.14  4 

3.9 

248  24 

.58  272 

7.4 

409  36 

.58  14 

2.2 

270  26 

.45  38 

5.1 

201  49 

.44  22 

7.7 

232  54 

.59  21 

7.2 

49 


TABLE  11(d) 

Mean  and  standard  deviation  of  slab  thickness  (km),  ratio  of  bottomside  TEC  to  topside  TEC, 
ionogram  count  and  foF2  (MHz)  for  Lattm  in  [30°  N,  6(f  N\,  Kp  >  2.5,  Rz  >  40. 


^  v  Local 

Month 

[0.S1 

[6.121 

[12.18] 

[18.24] 

January 

to 

March 

393  34 

.73  12 

2.8 

- 

231  46 

.48  19 

8.9 

- 

April 

to 

June 

298  108 

.49  18 

3.6 

- 

287  37 

.77  33 

8.5 

- 

July 

to 

September 

- 

306  35 

.40  27 

3.6 

251  26 

.47  22 

7.0 

- 

October 

to 

Oec ember 

- 

266  41 

.35  22 

8.5 

229  45 

.64  44 

7.8 

269  30 

.72  41 

5.2 

In  Table  12,  the  results  of  computing  slab  thickness  and  the  height  of 
the  base  of  the  F2  layer  have  been  compared  for  the  two  different  parabolas 
used  (with  and  without  ground  echoes) .  Of  the  198  determinations  made  (or 
ionograms  observed),  from  February  1965  and  1967,  a  systematic  selection 
(every  18th  determination)  has  been  displayed. 

The  empirical  model  of  slab  thickness,  Th,  displayed  in  COMPUTER  OUT¬ 
PUT  3  is  briefly  described.  The  result  of  sorting  the  10,000  ALOUETTE  NH 
profiles  into  384  cells  yielded  132  points  at  which  slab  thickness  is  known: 

^Lti’  Mi’  Lattmi»  KPi*  Lon81*  Th^,  i  =  1(1)132 

The  majority  of  these  points  occurred  in  the  two  most  northeast 
geographic  zones  since  that  was  where  the  majority  of  the  10,000  profiles 
occurred.  The  252  =  (384  -  132)  empty  cells  correspond  to  non-existent 
keyed  records  on  the  file  THICK,  (see  also  Appendix  A). 


These  points  were  then  fitted  to  a  model  of  the  form: 


51 


The  actual  6  independent  variables  used  were: 

x^  =  irLt/12  Lt  in  [0,24] 

x2  =  n (2M-1) /12  M  in  [.5,12.5] 

(integer  months  denote  mid-months) 

x^  =  Lattm/90  Lattm  in  [0,90] 

x^  =  Kp/9.3  Kp  in  [0,9] 

=  (Long  +  180) /140  Long  in  [-180°, -40° ] 

x,  =  Rz/120  Rz  in  [0,150] 

o 

The  statistically  significant  coefficients,  b.:  j=0(l)24,  were  computed  by 

the  subroutines  STEP1  and  STEP2.  In  this  case,  the  actual  coefficients  were: 

b  =  166.05  b,  =  -38.462  b.,  -  32.407 

o  1  3 

b5  =  63.112  b9  =  7.9556  b  =  -29.452 

b^  =  57.693  b._  =  100.47  b0/  =  119.29 

16  17  24 

and  the  root  mean  square  error  for  the  errors  in  predicted  thickness  was  34 
km.  As  a  matter  of  interest,  all  the  variables  except  longitude  appeared  in 
the  model  (3.18) . 

The  subroutine  SLABTH  in  the  main  prediction  program  PREDSYS  actually 
computes  slab  thickness  for  the  user-specified  conditions  and  for  each  hour 
of  local  time. 

In  conclusion,  the  following  remarks  pertaining  to  these  illustrations 
are  made: 

(a)  More  profile  data  is  needed  to  cover  the  entire  sunspot  range 
adequately.  These  data  exist  (see  ISIS  1  and  2,  Table  8)  where 
one  observes  that  there  are  approximately  70,000  more  profiles, 
in  the  North  American  sector,  which  may  be  used.  These,  combined 
with  the  early  Alouette  1  profiles  would  give  a  fairly  extensive 
data  base  of  approximately  130,000  ionograms  (=  50,000  +  10,000  + 
70,000).  While  this  is  still  only  about  one  third  of  the  size  of 
the  ALOSYN  data  base  used,  it  would  permit  a  finer  cell  division 
(than  used  in  Tables  10  -  11)  and  also  fill  in  some  of  the 
seasonal  gaps  observed  in  the  latter. 

(b)  Even  with  an  extended  data  base,  there  still  exists  the  uncertainty 
of  determining  the  profile  (and  TEC)  on  the  bottomside  of  the  F2 
layer,  from  topside  data  only.  While  it  may  not  be  feasible  to 
use  the  "ground  echo"  approach  everywhere  (recall  reasons  in 
section  3.2.2),  there  do  exist  extensive  bottomside  synoptic  data 
(scaled  and  recorded  on  magnetic  tape)  containing  ionospheric 
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characteristics  such  as  foF2,  foFl,  foE,  h'F2,  h'F,  M(3Q00)F2, 
and  M(3000)F1  for  ground  stations  located  at  Resolute,  Churchill, 
Winnipeg-Kenora,  Ottawa,  and  St.  John's,  at  least.  (See  Tapes 
5E08,  4B32,  and  4B33  in  Appendix  A).  These  data  would  certainly 
give  at  least  a  "representative"  parabola,  or  combination  of 
parabolas,  for  the  bottomside  profile,  using  the  classical  Abel 
equation  (see  [8]  of  Bibliography)  relating  scaled  ionograms  and 
electron  density  profiles. 


4.  MAIN  TROUGH  REFINEMENTS 
4.1  DESCRIPTION  OF  THE  PROCEDURE 

Beginning  with  Muldrew  [17],  extreme  troughs  have  been  consistently 
observed  in  the  spatial  electron  density,  particularly  in  the  sunset  to 
sunrise  period  of  the  day  and  at  mid  and  high  geomagnetic  latitudes.  In 
addition,  regular  movements  of  the  trough  have  been  observed  (equatorward 
with  increasing  magnetic  activity  and  local  time  tending  towards  midnight). 
Tulunay  and  Sayers  [28],  using  data  from  the  Ariel  III  satellite  launched  in 
1967,  compiled  a  statistical  summary  of  trough  crossings  over  an  11  month 
period.  In  particular,  they  noted  that  most  of  the  troughs  were  observed  in 
the  geomagnetic  latitude  range  [50°N,  70°N].  In  Pike  [22],  there  is  a  set 
of  linear  equations  for  predicting  foF2  from  a  "reference  latitude  coordinate 
system".  These  are  for  a  very  limited  period  of  time  (December  1971  ISIS  2 
data).  Halcrow  and  Nisbet  [10]  produced  a  trough  model  using  30,000  profiles 
of  Alouette  1  and  2  in  the  time  period  1962-70.  The  output  from  their  model 
was  a  multiplicative  factor,  depending  on  invariant  latitude,  Kp,  Lt,  and  Rz 
which  could  be  applied  to  the  Nmax  (or  foF2)  as  predicted  by  the  C.C.I.R. 
model,  [5],  Essentially  they  defined  it  as  unity  at  the  north  and  south 
edges  of  the  trough  and  0.25  at  the  centre,  for  Nmax. 

The  procedure  adopted  here  is  essentially  to  use  the  ALOSYN  model  to 
locate  the  trough  for  a  given  longitude,  month,  local  time,  sunspot  number, 
and  magnetic  activity.  The  changes  in  foF2  between  the  trough  edges  and  the 
centre  will  then  be  modified  to  agree  with  Halcrow  and  Nisbet's  observations 
to  partially  account  for  the  smoothing  that  has  occurred  due  to  the  averaging 
of  ionogram  parameters  in  the  cells. 

More  precisely,  the  procedure  is  as  follows: 

(1)  The  ALOSYN  model  is  used  to  generate  an  array  of  points  (Lattm^, 
foF2|c)  in  [45°,  75°  ]  (geomagnetic),  k  =  1(1)13,  for  a  given 
longitude,  local  time,  Kp,  month  and  Rz.  The  model  actually 
gives  pairs  at  every  5  degrees  of  latitude  -  additional  inter- 
polants  have  been  filled  in  using  a  cubic  spline  technique. 

(2)  Locate  the  trough  centre,  M;  i.e.  that  point  (Lattnty,  foF2M) 
where  foF2  takes  its  absolute  minimum  in  [45°,  75°]. 

(3)  The  north  and  south  edges,  N  and  S  respectively,  are  taken  to  be 
the  nearest  points  of  inflection  about  the  centre  M,  i.e.  (Lattm^, 
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foF2N)  and  (Lattms,  foF2g)  respectively.  This  is  the  reason  why 
interpolants  were  generated  above;  recall  that  the  trough  is  about 
5  degrees  wide  (north  to  south).  See  also  [10,17,19]. 

(4)  Finally,  the  corrected  foF2  at  the  centre  is  defined  by: 

f°F2c(M)  =  ™aX  ^  mln(foF2N»  f °F2g) ,  foF2%]  (4.1) 

where  foF2%  is  a  5-percentile  which  is  also  obtainable  from  the 
ALOSYN  model.  The  5-percentile  of  foF2  is  taken  from  the  distri¬ 
bution  of  foF2  at  the  point  (lattm^,  foF2M) .  An  example  of  an 
averaged  and  corrected  trough  is  now  presented. 


4.2  ILLUSTRATION 

In  Figures  13(a)-(d),  profiles  of  foF2  against  latitude  have  been 
displayed  for  the  indicated  conditions  (January  -  March,  1963  has  an  Rz  of 
approximately  29).  Unstandardized  (1963)  data  have  been  used  in  these  plots. 

All  the  plots  are  for  the  eastern-most  longitude  zone  [-90°,  -45°]. 

For  a  given  figure,  the  4  profiles  represent  the  4  local  time  zones  for  the 
month.  For  each  curve  the  local  time  can  range  up  to  5  to  6  hours  over  the 


Figure  13(a).  Illustration  of  the  main  trough  for  low  sunspot  number,  longitude  in  [9<fW,  45°  (V]  and  local  time 

in  Table  13,  January.  1963,  Kp  <  1.5. 


CO-GEOMAGNETIC  LATITUDE  (•) 
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entire  latitude  range.  For  convenience,  these  local  time  ranges  are  given  in 
Table  13.  Some  evidence  of  a  trough  (near  a  co-geomagnetic  latitude  of  30°) 
appears  in  all  the  figures  except  the  March  1963  midday  curves. 

In  Figure  13(e),  the  trough  correction  procedure  is  illustrated  for  the 
March  low  magnetic  activity  trough  which  appears  to  begin  to  form  at  about  19 
to  23  hours  local  time.  (One  should  observe  also  Figure  13(c)  and  note  that 
this  trough  has  become  more  accentuated  at  approximately  21  to  2  hours,  local 
time.  The  actual  step  by  step  trough  correction  is  now  illustrated.  In  the 
neighbourhood  of  the  trough  the  following  interpolants  have  been  generated: 


k 

Lattm 

f  oF2 

foF2% 

1 

75.0 

3.0 

2 

72.5 

3.0 

1.7 

3 

70.0 

2.9 

4 

67.5 

2.8 

1.7 

5 

65.0 

2.7 

6  N 

62.5 

2.6 

1.9 

7 

60.0 

2.3 

8  M 

57.5 

2.3 

1.1 

9 

55.0 

2.8 

10  S 

52.5 

3.5 

2.9 

11 

50.0 

3.8 

12 

47.5 

3.9 

2.4 

13 

45.0 

3.9 

as  well  as  a  5%  -  tile  of  foF2,  from  the  actual  cell  histograms.  The  pairs 
used  for  the  north  and  south  edge  as  well  as  the  centre  have  also  been 
appropriately  labelled.  Thus  the  corrected  foF2  at  the  centre  becomes, 
using  (4.1): 

foF2  =  max  [!j  min(2.6,  3.5),  1.1]  =  1.3  MHz 
The  centre  is  still  assumed  to  be  located  at  Lattm  =  57.5°N. 

In  conclusion,  the  prediction  system  (PREDSYS  and  associated  XEQ  file, 
i.e.  execute  file  in  Honeywell  CP-V  system,  EXPREDSYS)  has  been  set  up  so  as 
to  allow  the  user  to  optionally  make  trough  corrections,  by  hand,  using  the 
above  illustration  as  a  guide.  In  other  words,  after  one  has  made  sets  of 
predictions  for  several  latitudes  (at  most  19),  a  file  will  exist  containing 
grids  of  the  characteristics  (foF2,  spread,  and  slab  thickness  in  local  time 
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TABLE  13 

Local  time  code  for  Figures  13(a)  to  13(e) 


LOCAL  TIME 
CELL  INDICES 

GRAPHIC 

SYMBOL 

LOCAL  TIME  RANGE  (Hours) 

(ALL  LATITUDES) 

(LT1.LT2) 

JANUARY 

MARCH 

(1,1) 

(21,2| 

[13,18] 

(1,2) 

[5,10] 

[21,21 

(2,1) 

[19,24] 

[11,161 

(2,2) 

13,7] 

[19,23] 

and  latitude).  The  user  is  then  permitted  to  edit  this  file  (using  the  EDIT 
processor  of  the  Honeywell  Sigma  9  computer)  in  order  to  make  trough 
corrections,  fill  in  "missing  values"  resulting  from  the  lack  of  cells  in 
ALLMON,  and  even  correct  "obvious"  bad  predictions.  To  aid  in  this,  the  user 
may  access  those  cells  that  are  available  (see  COMPUTER  OUTPUT  4),  and  also 
access  files  of  preliminary  interpolates  (see  section  2.3.2),  which  are  also 
produced  by  the  prediction  system  at  the  user's  option. 

Upon  completion  of  this  editing  task,  the  user  may  resume  use  of  the 
XEQ  file  EXPREDSYS,  to  display  the  contour  maps. 
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APPENDIX  A 


Glossary  of  Files  and  Tapes 

In  this  Appendix,  the  files  associated  with  each  section  of  the  report 
are  listed  (in  application  order)  and  briefly  described.  Each  such  description 
is  prefixed  by  the  code: 

FILENAME,  TYPE,  ORG. 
where  (1)  FILENAME  is  the  name  of  file 

(2)  TYPE  is  one  of  the  letters  X,B,F,  or  D  corresponding  to  XEQ,  BATCH, 
FORTRAN  source,  or  data  files. 

(3)  ORG  is  one  of  K3,  Kll,  or  C  describing  the  key  structure  of  the  file. 
Here  K3  denotes  files  with  3  byte  (EDIT)  keys,  Kll  with  11  byte 
(direct  access)  keys,  and  C  for  sequential  files  (no  keys).. 

Any  F  file  is  compiled  on  the  file  BFILENAME  in  single  precision,  unless 
the  contrary  is  stated. 

Any  compiled  version  of  an  F  file  is  loaded  with  LYNX  and  the  file  BMISC, 
unless  the  contrary  is  stated. 

In  some  cases,  for  a  related  set  of  a  files,  two  numbers  denoted  by  XX 
will  be  suffixed  to  the  file  name  and  will  be  defined  at  that  point. 

The  file,  EXFILENAME,  will  denote  an  X  file  for  the  associated  F  file, 
FILENAME.  The  source  listings  of  all  X,B  and  F  files  are  available  on  request. 
See  Tape  5E47  below. 

The  first  file  is  of  general  purpose  and  applies  to  the  entire  report. 
MISC,F,K3. 

All  miscellaneous,  non-standard  subroutines  not  displayed  explicitly 
elsewhere  are  here.  It  should  be  compiled  with  the  batch  compile  (BC)  option. 


Section  2.2.1 

CRALOMON , X , K3 . 

This  file  controls  the  transfer  of  the  original  ALOSYN  tapes  to  disk. 


DXX , D , K3 .  XX  -  62(1)68 

These  7  files  are  built  by  the  user  in  CRALOMON.  Each  contains  1  record 
with  2  integers  which  are  used  for  positioning  an  ALOSYN  tape  to  the  month  to 
be  processed.  See  also  IONMODELS. 


DM0N.D.K3. 

This  file  is  built  in  CRALOMON  and  consists  of  1  record  with  8  integers 
to  be  used  for  defining  the  cell  alignment.  See  also  IONMODELA. 

ALOMONXX , D , C .  XX  =  62(1)68 

These  7  files  each  consist  of  the  same  month  of  ALOSYN  data,  for  each 
year  XX. 

CXX,D,C.  XX  =  62(1)68 

Each  of  these  files  contain  5  records  with  certain  quality  statistics 
concerning  the  respective  ALOMONXX.  See  also  IONMODELS. 

IONMODELS , F , K3 . 

This  program  produces  ALOMONXX  and  CXX,  one  year  at  a  time. 

DUMFILE,-,K3. 

Records  9  to  28  are  PCL  commands  used  for  removing  and  accessing  ALOSYN 
tapes  as  they  are  processed. 

CNSTRPCL,F,K3. 

This  is  used  for  copying  the  last  3  records  of  ALOMONXX.  The  actual 
number  of  records  of  ALOMONXX  is  given  by  the  first  record  of  the  appropriate 
CXX.  It  should  be  compiled  with  the  BC  option. 

AXX,B,K3.  XX  =  62(1)68 

These  7  files  control  the  production  of  DUMXX,  respectively.  They  are 
"batched"  by  CRALOMON. 

BXX,B,K3 .  XX  =  62(1)68 

These  7  files  control  the  production  of  MONXX,  respectively.  They  are 
"batched"  by  CRALOMON. 

Section  2.2.2 

DUMXX, D,C.  XX  =  62(1)68 

These  7  files  contain  the  cell  indices,  input  averages,  and  extreme  foF2 
values  for  the  ionograms  of  ALOMONXX. 
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IONMODELA , F , K  3 

This  produces  DUMXX  from  ALOMONXX,  one  year  at  a  time. 

MONXX ,  D ,  C .  XX  =  62(1)68 

These  7  files  are  the  respective  DUMXX  updated  with  the  histograms  of 
foF2  and  spread  F. 

IONMODELB , F , K3 . 

This  produces  MONXX  given  an  ALOMONXX  and  a  DUMXX. 


REDUALO , F , K3 . 

The  orbital  position  parameters  of  the  satellite  are  corrected  on  ALOMONXX 
using  this  program,  (see  February,  1965). 


IONGACC , F , K3 . 

Individual  ionograms  from  user  specified  cells  are  accessed  from  ALOMONXX. 


SATTRACK,F,K3. 

This  creates  a  file  of  orbital  position  parameters  of  the  satellite,  for 
a  given  time  interval.  These  data  are  used  in  the  construction  of  Figure  4. 
See  also  EXSATTRACK. 


Section  2.2.3 
MERMON,X,K3. 

This  controls  the  merging  and  blocking  of  the  MONXX  files,  as  well  as 
the  production  of  ALLMON. 


DUMFILE,- ,K3 . 


Records  2  to  8  contain  merging  and  blocking  specification  parameters 
required  by  the  system  MERGE  processor.  Record  1  contains  two  integers  which 
are  required  as  input  for  IONMODELC. 


MON , D , C 


The  merged  version  of  all  the  MONXX  files 
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MONB96 ,D,C . 

The  blocked  version  of  MON.  Nyear  years  for  a  given  cell  are  blocked 
together  on  a  record  (in  increasing  year  order).  Each  record  has  96*Nyear 
bytes.  All  empty  cells  are  included,  in  order  to  facilitate  the  accession  of 
user-specified  cells.  These  files,  all  months,  have  been  stored  on  magnetic 
tape.  See  tape  listing  below. 


ZR6278,D,K3. 

This  consists  of  34  records;  the  odd  records  are  the  monthly  Zurich 
sunspot  numbers  for  each  year  of  1962-78,  while  the  even  records  are  the 
smoothed  monthly  Zurich  sunspot  numbers . 


YEAR,D,K3. 

This  contains  one  record  with  two  integers  specifying  the  month  and  the 
number  of  years  to  be  processed  by  IONMODELC.  See  also  record  1  of  DUMFILE. 


IONMODELC , F , K3 . 

This  standardizes  the  records  of  MONB96  to  the  average  sunspot  number 
of  the  cell. 


MONSTD.D.KU. 

A  keyed  file  corresponding  to  the  standardized  month  represented  by 
MONB96  and  produced  by  IONMODELC.  Only  those  cells  to  be  used  in  interpolation 
are  kept. 


ALLM0N,D,K11. 

The  concatenated  version  of  MONSTD  for  all  months.  The  maximum  allowable 
key  is  FFFFFF^g  or  16,777,215^0.  The  actual  maximum,  for  the  cell  arrangement 
described  in  this  report  is  10,368,000  («  12  x  864  x  1000).  Keys  less  than 

10,000,000  can  be  directly  displayed  using  the  EDIT  processor  (system).  Other¬ 
wise  the  system  processor,  KFM  (see  Communications  Research  Centre  User's  (COMUS) 
Manual  5-2.13)  must  be  invoked.  Note  that  the  keys  as  defined  in  FORTRAN  are 
in  the  range  [1,10368]. 


KEYKFM,F,K3. 

This  is  used  for  accessing  keys  of  ALLM0N  which  are  not  compatible  to 
EDIT.  See  also  EXKEYKFM . 


ACCMONB96,F,K3. 

This  is  used  for  accessing  cells  from  the  unstandardized  months,  M0NB96, 
upon  specification  of  the  cell  indices.  Statistics  of  the  characteristics  are 
also  computed. 

CELLCOUNT , F , K3 . 

This  is  for  examining  the  quality  of  the  cells  of  MONB96,  as  in  Table  5. 


CELLKEY.F.K3. 

This  is  used  for  the  construction  of  Table  6. 


CHRZGRAD , F , K3 . 

This  is  applied  to  MONSTD  in  order  to  smooth  foF2-Rz  gradients.  (See 
also  Figures  6(a)  and  6(b)). 


Section  2.2 

PREDSYS,F,K3. 

This  is  the  prediction  program  which  uses  keyed  records  of  ALLMON  to 
interpolate  and  predict  ionospheric  characteristics.  It  should  be  compiled 
with  the  BC  option,  as  it  contains  the  subroutines  REGRESS,  SLABTH,  and  DIVTEST. 
Also,  to  load,  LYNX  with  BSSTEP1  and  BSSTEP2. 


REGRESS, F,K3. 

This  is  used  to  prepare  the  arrays  for  the  subroutines  STEP1  and 

STEP2. 


SLABTH, F,K3. 

Slab  thickness  is  computed  for  the  user's  input  conditions,  using  the 
model  illustrated  in  equation  (3.18). 


DIVTEST, F,K3. 

Any  positive  integer  is  tested  for  divisibility  by  a  given  integer. 

This  facilitates  the  accession  of  the  6  half  months  of  data  with  2  periods  of 
local  time  each  which  are  used  in  the  model  (2.13). 
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PRDATA,D,K3. 

Each  record  contains  a  set  of  Input  conditions  for  making  predictions  at 
a  given  longitude,  latitude,  sunspot  number,  month,  and  magnetic  activity. 

See  also  tape  R033  for  an  example. 


ADJRZ ,D,C . 

Preliminary  interpolants  (characteristics  adjusted  to  the  user's  sunspot 
number)  are  displayed. 


ADJKP , D , C . 

The  (at  most)  12  triples  of  month  and  local  time  data  are  displayed. 


REGFIT,D,C. 

The  details  of  the  stepwise  regression  procedure  for  the  model  (2.13) 
are  displayed. 


SSTEP1  and  SSTEP2,F,K3 

These  contain  the  subroutines  STEP1  and  STEP2  which  perform  the  stepwise 
regression  for  fitting  a  characteristic  as  a  function  of  local  time  and  month. 


NEWPL0T,F,K3 

Contour  maps  of  characteristics  in  local  time  and  latitude  are  prepared. 
Compile  with  the  BC  option  to  include  the  subroutine  CPLOTOPT. 


CPLOTOPT , F , K3 . 

A  modified  version  of  the  MLIB's  contour  plotting  routine. 


NEWPL0T1 ,D,C . 

This  contains  the  contour  grids  produced  by  PREDSYS.  They  may  be 
(optionally)  edited,  using  EDIT,  before  being  input  to  NEWPLOT. 


EXPREDSYS,X,K3. 


This  controls  the  entire  procedure  described  in  section  2.3 
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Section  3. 2. 1 

NHA.D.C. 

A  file  consisting  of  original  N(h)  profiles  in  A  Format  (see  Table  8). 


NHB.D.C. 

A  file  consisting  of  original  N(h)  profiles  in  B  Format  (see  Table  8). 


NHAMES.D.C. 


A  file  consisting  of  original  N(h)  profiles  in  Ames  Format  (see  Table  8). 
KPFILE.D.C. 

A  file  such  that  a  record  consists  of  the  8  3-hourly  Kp  indices  together 
with  the  associated  year  and  day. 


NHXX,D,C . 

A  file  consisting  of  reduced  profiles  in  the  format  of  ALOUETTE  NH.  XX 
are  two  integers  denoting  the  satellite  number  in  [1,2]  and  the  volume  number 
in  [1,9]  (see  Table  8). 


REDNHA.F.K3. 


This  produces  a  particular  NHXX  from  a  file  NHA. 


REDNHB, F,K3. 

This  produces  a  particular  NHXX  from  a  file  NHB. 


AMESPR,F,K3. 

This  produces  a  particular  NHXX  from  a  file  NHAMES.  It  should  be 
compiled  with  the  "Batch  Compile”  (BC  option)  as  it  includes  the  subroutine 
FIELDAL. 


NHCSB.D.C. 

This  consists  of  all  the  NHXX;  concatenated,  sorted  (by  universal  time) 
and  blocked  (at  the  rate  of  10  profiles  per  record).  An  extra  byte  (garbage) 
has  been  added  to  each  profile  to  facilitate  further  buffering.  The  actual 
procedure  is  listed  here  in  the  execute  file  EXNHCSB. 
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EXNHCSB 
!PCL  DUMMY 1 

•SET  F : SORTOUT /NHCSB ;OUT 
! SET  F : SORTIN/NHC ; IN 
.•SET  M:SI/DUMMY2;IN 
•SORT 

where  DUMMY1  consists  of: 

C  NHXX , NHXX , . . . , NHXX  INTO  NHC 
END 

and  DUMMY2  consists  of: 

•REC(103,104) 

.KEYS (9, 15) 

. BLOCK ( 1 ,  10, DROP) 

Section  3.2.2 

THICKEX,F,K3. 

Profiles  from  ALOUETTE  NH  are  buffered  in  and  a  slab  thickness  is  defined 
for  each.  These  results  are  then  sorted  into  cells,  with  an  alignment 
specified  by  the  user,  and  displayed  in  the  form  of  a  direct  access  file. 


THICK, D,K11 

The  output  of  THICKEX. 


THICKPR,F,K3. 

THICK  is  read  in  a  sequential  fashion,  and  a  model  as  in  (3.18)  is 
constructed.  LYNX  with  BSSTEP1  and  BSSTEP2. 


Section  3.2.3 
NHCS,D,C. 

A  selection  of  (unblocked)  103  byte  profiles  of  ALOUETTE  NH. 
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ALO,D,C. 


A  portion  of  a  file  ALOMONXX  which  contains  the  ionograms  listed  in  NHCS. 
ALOSNHMAT , F , K3 . 

A  program  which  extends  the  records  of  NHCS  (to  132  bytes)  by  attaching 
the  foF2  and  related  parameters  of  ALOSYN.  Table  9  was  also  produced  by  this 
program. 


NHF0F2,D,C. 

This  is  the  output  file  produced  by  ALOSNHMAT. 


GRECHO , D , C . 

A  file  consisting  of  the  output  on  unit  108  of  RANDRUN.284201E  (updated 
May  3,  1974  [J.H.  Whitteker,  private  communications,  1978],  The  card  input 
(units  105  and  306)  essentially  consist  of  the  scaled  ground  echoes  (X  wave) 
of  the  ionograms  in  NHF0F2.  The  Number  density  contour  plots  on  unit  108  are 
not  included. 


GRECHOM.D.C. 

A  data  file  consisting  of  m(n+l)  records  where  m  is  the  number  of 
ionograms  of  NHFOF2  scaled  for  ground  echoes  and  n  is  the  number  of  pairs  of 
echoes  actually  scaled  for  each  ionogram.  The  extra  record  is  an  identifi¬ 
cation  record  (header)  for  each  ionogram. 


CGRECHO , F , K3 . 

The  program  which  produces  GRECHOM  from  the  files  NHF0F2  and  GRECHO. 


EAREC.F.K3. 

The  complete  profiles,  corresponding  to  those  of  NHCS,  are  input  from 
NHA  and  NHB  together  with  the  file  GRECHOM  to  compute  a  slab  thickness.  It 
should  be  compiled  with  the  BC  option  (all  non-standard  subroutines  are 
included).  FIELDG  of  ORBIT. 284201C  and  SOLVLE  of  .ML IB  are  required  if  a 
double  precision  version  is  desired.  The  input  files  GRECHOM,  NHB,  and  NHA 
used  for  Table  12  are  listed  on  tape  6E19,  see  below. 


Section  4 

TR0CRT,F,K3. 

A  profile  of  a  characteristic  against  latitude  is  created  (as  a  smooth 
curve,  using  the  subroutine  CUBSPL,  which  is  included)  from  a  grid  of  the 


characteristic,  in  the  format  of  Table  6.  Thus  TROCRT  should  be  compiled  with 
the  BC  option.  CUBSPL  also  uses  SOLVLE  which  may  be  found  in  the  source  file 
SSOLVLE,  for  single  precision  use. 


SS0LVLE,F,K3 

A  single  precision  version  of  the  program  SOLVLE,  (in  the  library 
programs  .MLIB  of  the  Honeywell  Sigma  9  computer)  for  solving  a  system  of 
linear  algebraic  equations. 

The  actual  tape  serial  numbers  used  in  this  report  are  now  listed. 

Unless  the  contrary  is  stated  these  tapes  are  all  9  track,  800  bytes  per  inch, 
and  free  tapes.  For  labelled  tapes,  the  serial  number  is  the  "Internal" 
serial  number,  i.e.  that  written  on  the  tape. 

1 .  ALOSYN  TAPES 


(a) 

6C47 

September  -  December, 

1962 

(b) 

6G03 

January 

-  April,  1963 

(c) 

4G48 

May 

-  December, 

1963 

(d) 

4G49 

January 

-  December, 

1964 

(e) 

4E26 

January 

-  December , 

1965 

(f) 

1G47 

January 

-  December, 

1966 

(g) 

1G48 

January 

-  December, 

1967 

(h) 

1G49 

January 

-  December, 

1968 

(i) 

3 BO 2  or  6E21 

MONB96, 

for  all  months. 

ALOUETTE  N(h)  TAPES 

(a) 

R541 

Alouette 

1,  Volume  1 

(b) 

R648 

Alouette 

1,  Volume  2 

(c) 

R034 

Alouette 

1,  Volume  3 

(d) 

R535 

Alouette 

1,  Volume  4 

(e) 

R568 

Alouette 

1,  Volume  5 

Alouette  1,  Volume  6 

Alouette  1,  Volume  7,  See  3(c),  EXNHAL1. 


(f)  R571 

(g)  R033 
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(h) 

R037 

Alouette 

2, 

Volume  1 

(i) 

R752 

Alouette 

2, 

Volume  2 

(j) 

R582 

Alouette 

2, 

Volume  3 

(k) 

R740 

Alouette 

2, 

Volume  4 

(1) 

R570 

Alouette 

2, 

Volume  5 

(m) 

R022 

Alouette 

2, 

Ames  1 

(n) 

R720 

Alouette 

2, 

Ames  2 

(o) 

R739 

Alouette 

2, 

Ames  3 

(P) 

R749 

Alouette 

2, 

Ames  4 

(q) 

5E46  or  3E30 

ALOUETTE 

NH 

(or  NHCSB) 

(r) 

6E19 

GRECH0M, 

(2427  records) 

NHB 

NHA 


,  (1201  records) 
,  (  501  records) 


3.  MISCELLANEOUS  TAPES 

(a)  3B04 

(b)  5E47 


Kp  ,  January  1962  -  December  1970. 


Labelled  tape  (account  284201R, GLADWIN) 

This  contains  all  the  files  which  are  available 
on  request. 


(c)  R033 


Labelled  tape  (account  284201R, GLADWIN) 
This  contains  the  data  files 
ALLMON,  CONTOUR,  EXNHAL1,  FQJAN63, 
FOFEB63,  FOMAR63,  PRDATA,  THICK,  ZR6278. 


(d)  5E08 

(e)  4B32 


(f)  4B33 

(g)  R603 


(h)  R579 


Bottomside  synoptics:  January,  1969  -  December,  1971 

Bottomside  synoptics:  January,  1972  -  February,  1973 

Bottomside  synoptics:  January,  1973  -  December,  1973 

Alouette  1  (103  byte)  Profiles  ;  NHCS 
Alouette  1  (103  byte)  Profiles  ;  NHCS 
Alouette  2  (103  byte)  Profiles  ;  NHCS 
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(i)  R565 


APPENDIX  B 


List  of  Notations 

Individual  numbers  are  labelled  by  small  English  letters,  possibly 
subscripted  in  the  case  of  a  serially  related  set  of  numbers. 

File  names  (FORTRAN  programs  or  data  files)  are  labelled  with  capital 
letters. 

The  following  notations  are  also  used  for  denoting  sets  of  numbers  or 
files. 

(a)  A  -  [a,b]  :  The  set  of  numbers  lying  between  a  and  b  Including 
the  end  points,  with  the  additional  convention  that  a  round 
bracket  denotes  exclusion  of  the  end  point,  (a  <  b) . 

(b)  A  =  Ja,b,c,...j  :  The  set  of  numbers  consisting  of  a,  b,  c  ...  . 

(c)  If  the  elements  of  A  are  serially  related  the  following  is  used: 

A  -  |xt  =  i  -  n(r)n[  ...  .  x^J 

(d)  The  following  are  for  combination  of  sets: 

AUB:  The  set  consisting  of  the  elements  in  the  set  A  or  the 
set  B. 

AflB:  The  set  consisting  of  the  elements  in  the  set  A  and  the 
set  B. 

AxB:  The  set  of  all  combinations  (pairs)  of  elements;  the  first 
from  set  A  and  the  second  from  set  B. 

A-B:  The  set  of  elements  in  the  set  A  but  not  in  B. 


Chapters,  sections,  and  subsections  are  ordered  lexicographically  by 
a  triple  of  integers.  Equations  are  numbered  consecutively  within  a 
given  chapter,  by  a  pair  of  numbers. 
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NOTES,  ADDENDA,  AFTERTHOUGHTS 
on  CRC  Report  No.  1337 
C.J.  Gladwin 
January  27  1981 


,IJ  ■ 

rru| 

i  .  vwtl 


1.  Page  7:  There  "may  be"  more  data  for  October,  1964.  Note  in  Table  2 

that  the  number  of  blocked  records  for  this  month  is  considerably 
less  than  October  1963  and  October  1965.  See  also  a  tape  among 
the  original  ALQSYN  tapes  "missing  data  for  October  1964". 

However,  this  may  already  have  been  included  in  the  main  1964 
ALQSYN  tape.  In  any  case  to  do  a  given  month,  one  would  apply  the 
XEQ  and  BATCH  control  files  CRALOMON,  AXX,  BXX,  MERMON  to  an 
appropriate  set  of  ALQSYN  records,  AL00CT64.  (i.e.  AXX  *s  A64,  BXX 

=  B64,  ALOMONXX  =  ALOOCT64) 


2.  Page  16:  The  "format"  of  MONB96  *  SEPB96  (See  also  tape  3B02,  for 
1  example)  does  not  appear  exactly  as  in  COMPUTER  OUTPUT  1. 

Actually,  every  7  lines  (years)  are  written  as  a  blocked  record 
with  7  x  96  *  672  bytes.  (The  order  and  spacing  of  the  various 
numbers,  within  a  line,  is  preserved  though). 


3<  Page  24:  The  (median)  foF2  can  be  calculated  in  two  "equivalent"  ways, 
given  the  assumptions  for  the  distribution.  (See  also  Figure  7) • 


(a)  As  in  equation  (2.10),  Page  24. 
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(b)  The  median  of  the  standardized  foF2  distribution  can  be 

fouid  and  then  that  number  can  be  adjusted  to  the  user's  sunspot  nunbe 


f\  -f  ^ 


2 


S  j,. 


Method  (a)  has  been  used  in  the  main  predictions 
program.  To  implement  method  (b),  the  following  adjustment 
would  have  to  be  made  (in  PREDSYS): 

(i)  Delete  line  161,  "F0F2  =  A«USS+B" 

(ii)  Insert  lines  161,  161.1 

CALL  XFRACT(SFDISTR, 10,NGRAND,SMA,BIG, .5.F0F2) 

F0F2  =  F  OF  2+A* ( USS-SSA V ) 

The  "advantage"  of  this  second  approach  is  that  the 
median  and  percentiles  of  foF2  will  more  likely  be  correctly  ordered: 

i.e.  foF2  <  foF2  for  any  0  <a.  <a2  <  i 
O]  a2  1 

where  foF2  is  a  particular  percentile  or  median. 

These  characteristics,  however,  will  still  not  necessarily  be 
correctly  ordered  with  100%  certainty,  since  separate  (independent) 
models  of  the  form  (2.13)  are  fitted  to  each  characteristic. 

The  "disadvantage"  of  this  second  approach  is  that  one  is 
masking  cases  where  the  assumptions  of  Figure  7  are  not  too  good. 


4.  Page  26:  In  (2.12), 
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5.  P^e  51:  Delete  degrees  ("°")  in  x5  below. 

6,  Page  57:  Some  additional  hints  on  editting  the  output  of  PREDSYS,  namely 

NEWPL0T1. 

(a)  Examine  the  desired  portions  of  ADJRZ  to  see  if  the  order  of 
Lagrange  interpolations  of  foF2  with  Kp,  for  example,  agrees  with 
the  appropriate  (divided)  differences. 

A 

e.g.  Long  «  70PW,  Lattm  »  57°N,  Lt»  2200  hrs 

January,  Rz  =  30  (See  also  COMPUTER  OUTPUT  4) 


Kp 

fpF2  »  aRz  +  b 

AfoF2/6Kp 

0.52 

1.87 

0.13 

2.47 

2.12 

-0.17 

4.32 

1.80 

The  fourth  Kp  cell  only  had  1  year  of  data  (See  also  COMPUTER  OUTPUT  4), 
so  no  coefficients  a  or  b  were  available.  Thus,  Lagrange  interpolation 
of  order  2  (i.e.  a  parabola)  was  used.  In  this  case  second  differences 
"appear  constant"  as  they  should. 


I 


(b)  The  user  nay  wish  to  ignore  cells  of  ALLMON  which  yield  large 

i  ■  , 

sunspot  slopes,  "a".  (There  are  126  cells  with  a  >  0.1,  see  •Iso 
remarks  on  page  17). 

Also  he  could  agree  to  redefine  any  "a"  larger  than  0.1,  for 
example,  as  0.1.  This  latter  approach  is  probably  much  more  dangerous 
since  the  large  "a"  may  only  be  the  result  of  computing  with  2  years'  data 
at  essentially  the  sane  sunspot  number. 

Appropriate  restrictions  may  be  inserted  in  PRED5YS  between  lines 
156  and  157. 

• 

(c)  No  attempt  was  made  in  the  progrjam,  PREDSYS ,  to  mask  predictions 
which  yield  negative  values,  (for  instance  foF2  percentiles  for  low 
sunspot  nunbers,  or  troughs).  It  is  left  to  the  user's  "intuition" 
to  replace  such  values. 

« 

,  .  /  (For  exanple,  one  could  simply  use  "the  smallest  positive  value 

at  the  nearest  local  time".) 

7«  Page  74*75:  The  only  tape  serial  numbers,  with  any  real  meaning  to  the 
system,  are  the  nunbers  of  the  2  labelled  tapes. 


